Coherent searches for periodic gravitational waves from unknown isolated sources and 
Scorpius X-1: results from the second LIGO science run. 
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We carry out two searches for periodic gravitational waves using the most sensitive few hours of 
data from the second LIGO science run. Both searches exploit fully coherent matched filtering and 
cover wide areas of parameter space, an innovation over previous analyses which requires considerable 
algorithm development and computational power. The first search is targeted at isolated, previously 
unknown neutron stars, covers the entire sky in the frequency band 160-728.8 Hz, and assumes a 
frequency derivative of less than 4 x 10"^" Hz/s. The second search targets the accreting neutron 
star in the low-mass X-ray binary Scorpius X-1 and covers the frequency bands 464-484 Hz and 604- 
624 Hz as well as the two relevant binary orbit parameters. Due to the high computational cost of 
these searches we limit the analyses to the most sensitive 10 hours and 6 hours of data respectively. 
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Both searches look for coincidences between the Livingston and Hanford 4-km interferometers. 
Given the limited sensitivity and duration of the analyzed data set, we do not attempt deep follow- 
up studies. Rather we concentrate on demonstrating the data analysis method on a real data 
set and present our results as upper limits over large volumes of the parameter space. For isolated 
neutron stars our 95% confidence upper limits on the gravitational wave strain amplitude range from 
6.6 X 10"^^ to 1 X 10"^^ across the frequency band; For Scorpius X-1 they range from 1.7 X 10"^^ to 
1.3 X lO"'^^ across the two 20-Hz frequency bands. The upper limits presented in this paper are the 
first broad-band wide parameter space upper limits on periodic gravitational waves from coherent 
search techniques. The methods developed here lay the foundations for upcoming hierarchical 
searches of more sensitive data which may detect astrophysical signals. 

PACS numbers: 04.80.Nn, 95.55.Ym, 97.60.Gb, 07.05. Kf 



I. INTRODUCTION 



Rapidly rotating neutron stars are the most likely 
sources of persistent gravitational radiation in the fre- 
quency band « 100 Hz — 1 kHz. These objects may 
generate continuous gravitational waves (GW) through 
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a variety of mechanisms, including nonaxisymmetric dis- 
tortions of the star 13, la, ^0,13, velocity perturbations 
in the star's fluid Q^laQj and free precession |3, 0]. 
Regardless of the specific mechanism, the emitted signal 
is a quasi-periodic wave whose frequency changes slowly 
during the observation time due to energy loss through 
gravitational wave emission, and possibly other mecha- 
nisms. At an Earth-based detector the signal exhibits 
amplitude and phase modulations due to the motion of 
the Earth with respect to the source. The intrinsic grav- 
itational wave amplitude is likely to be several orders 
of magnitude smaller than the typical root-mean-square 
value of the detector noise, hence detection can only be 
achieved by means of long integration times, of the order 
of weeks to months. 

Deep, wide parameter space searches for continuous 
gravitational wave signals are computationally bound. 
At fixed computational resources the optimal sensitiv- 
ity is achieved through hierarchical search schemes [Tol 
Ull lia | . Such schemes alternate incoherent and coherent 
search stages in order to first efficiently identify statisti- 
cally significant candidates and then follow them up with 
more sensitive, albeit computationally intensive, meth- 
ods. Hierarchical search schemes have been investigated 
only theoretically, under the simplified assumption of 
Gaussian and stationary instrumental noise; the compu- 
tational costs have been estimated only on the basis of 
counts of floating point operations necessary to evalu- 
ate the relevant detection statistic and have not taken 
into account additional costs coming e.g. from data in- 
put/output; computational savings obtainable through 
efficient dedicated numerical implementations have also 
been neglected. Furthermore, general theoretical investi- 
gations have not relied on the optimizations that can be 
introduced on the basis of the specific area in parameter 
space at which a search is aimed. 

In this paper we demonstrate and characterize the co- 
herent stage of a hierarchical pipeline by carrying out 
two large parameter space coherent searches on data col- 
lected by LIGO during the second science run with the 
Livingston and Hanford 4-km interferometers. The sec- 
ond LIGO science run took place over the period 14 Feb. 
2003 to 14 Apr. 2003. As we will show, this analysis 
requires careful tuning of a variety of search parameters 
and implementation choices, such as the tilings of the pa- 
rameter space, the selection of the data, and the choice of 
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the coincidence windows, that arc difficult to determine 
on purely theoretical grounds. This paper complements 
the study presented in |0| where we reported results ob- 
tained by applying an incoherent analysis method to 
data taken during the same science run. Furthermore, 
here we place upper limits on regions of the parameter 
space that have never been explored before. 

The search described in this paper has been the 
test-bench for the core science analysis that the Ein- 
stein@home fl^ project is carrying out now. The de- 
velopment of analysis techniques such as the one de- 
scribed here, together with the computing power of 
EinsteinQhome in the context of a hierarchical search 
scheme, will allow the deepest searches for continuous 
gravitational waves. 

In this paper the same basic pipeline is applied to and 
tuned for two different searches: (i) for signals from iso- 
lated sources over the whole sky and the frequency band 
160 Hz - 728.8 Hz, and (ii) for a signal from the low-mass 
X-ray binary Scorpius X-1 (Sco X-1) over orbital param- 
eters and in the frequency bands 464 Hz ~ 484 Hz and 
604 Hz - 624 Hz. It is the first time that a coherent anal- 
ysis is carried out over such a wide frequency band, using 
data in coincidence and (in one case) for a rotating neu- 
tron star in a binary system; the only other example of 
a somewhat similar analysis is an all-sky search over two 
days of data from the Explorer resonant detector over a 

0. 76 Hz band around 922 Hz [lllllll^. 

The main scope of the paper is to illustrate an analysis 
method by applying it to two different wide parameter 
spaces. In fact, based on the typical noise performance 
of the detectors during the run, which is shown in Fig.^ 
and the amount of data that we were able to process in 
w 1 month with our computational resources (totalling 
about 800 CPUs over several Beowulf clusters) we do 
not expect to detect gravitational waves. For isolated 
neutron stars we estimate (see Section IIIII for details) 
that statistically the strongest signal that we expect from 
an isolated source is 4 x 10^^'* which is a factor ^ 20 
smaller than the dimmest signal that we would have been 
able to observe with the present search. For Scorpius X- 

1, the signal is expected to have a strength of at most ~ 
3 X 10~^^ and our search is a factor ~ 5000 less sensitive. 
The results of the analyses confirm these expectations 
and we report upper limits for both searches. 

The paper is organized as follows: in Section ^ we 
describe the instrument configuration during the second 
science run and the details of the data taking. In Sec- 
tion ^H] we review the current astrophysical understand- 
ing of neutron stars as gravitational wave sources, in- 
cluding a somewhat novel statistical argument that the 
strength of the strongest such signal that we can expect 
to receive does not exceed /ig""' « 4 x 10~^^. We also 
detail and motivate the choice of parameter spaces ex- 
plored in this paper. In Section llVI we review the signal 
model and discuss the search area considered here. In 
Section^we describe the analysis pipeline. In Section lVll 
we present and discuss the results of the analyses. In Sec- 



tion we recapitulate the most relevant results in the 
wider context and provide pointers for future work. 



II. INSTRUMENTS AND THE SECOND 
SCIENCE RUN 

Three detectors at two independent sites comprise the 
Laser Interferometer Gravitational Wave Observatory, or 
LIGO. Detector commissioning has progressed since the 
fall of 1999, interleaved with periods in which the ob- 
servatory ran nearly continuously for weeks or months, 
the so-called "science runs". The first science run (SI) 
was made in concert with the gravitational wave detec- 
tor GEO600; results from the analysis of those data were 
presented in [Tsl l2l[ |2^ |2^ , while the instrument sta- 
tus was detailed in |2J| . Significant improvements in the 
strain sensitivity of the LIGO interferometers (an order 
of magnitude over a broad band) culminated in the sec- 
ond science run (S2), which took place from February 14 
to April 14, 2003. Details of the S2 run, including de- 
tector improvements between SI and S2 can be found in 
O, Section IV of HI, and Section II of and |2|. 

Each LIGO detector is a recycled Michelson interfer- 
ometer with Fabry-Perot arms, whose lengths are defined 
by suspended mirrors that double as test masses. Two 
detectors reside in the same vacuum in Hanford, WA, one 
(denoted HI) with 4-km armlength and one with 2-km 
(H2), while a single 4-km counterpart (LI) exists in Liv- 
ingston Parish, LA. Differential motions are sensed inter- 
fcrometrically, and the resultant sensitivity is broadband 
(40 Hz - 7 kHz), with spectral disturbances such as 60 Hz 
power line harmonics evident in the noise spectrum (see 
Fig.^l. Optical resonance, or "lock", in a given detector 
is maintained by servo loops; lock may be interrupted by, 
for example, seismic transients or poorly conditioned ser- 
vos. S2 duty cycles, accounting for periods in which lock 
was broken and/or detectors were known to be function- 
ing not at the required level, were 74% for HI, 58% for 
H2, and 37% for LI. The two analyses described in this 
paper used a small subset of the data from the two most 
sensitive instruments during S2, LI and HI; the choice 
of the segments considered for the analysis is detailed in 
Sec. ED 

The strain signal at the interferometer output is recon- 
structed from the error signal of the feedback loop which 
is used to control the differential length of the arms of 
the instrument. Such a process — known as calibration — 
involves the injection of continuous, constant amplitude 
sinusoidal excitations into the end test mass control sys- 
tems, which are then monitored at the measurement error 
point. The calibration process introduces uncertainties in 
the amplitude of the recorded signal that were estimated 
to be 11% during S2 |23|. In addition, during the run 
artificial pulsar-like signals were injected into the data 
stream by physically moving the mirrors of the Fabry- 
Perot cavity. Such "hardware injections" were used to 
validate the data analysis pipeline and details are pre- 
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FIG. 1: Typical one-sided amplitude spectral densities of de- 
tector noise during the second science run, for the three LIGO 
instruments. The solid black line is the design sensitivity for 
the two 4-km instruments LI and HI. 

sented in Sec. IV HI 

III. ASTROPHYSICAL SOURCES 

We review the physical mechanisms of periodic gravi- 
tational wave emission and the target populations of the 
two searches described in this paper. We also compare 
the sensitivity of these searches to likely source strengths. 

A. Emission mechanisms 

In the LIGO frequency band there are three predicted 
mechanisms for producing periodic gravitational waves, 
all of which involve neutron stars or similar compact ob- 
jects: (1) nonaxisymmetric distortions of the solid part 
of the star 

QJ,I1MJ1L(2) unstable r-modes in the fluid 
part of the star IS Qli ^^^^ (3) free precession of the 
whole star 

We begin with nonaxisymmetric distortions. These 
could not exist in a perfect fluid star, but in realistic 
neutron stars such distortions could be supported either 
by elastic stresses or by magnetic fields. The deformation 
is often expressed in terms of the ellipticity 



which is (up to a numerical factor of order unity) the 
771 = 2 quadrupole moment divided by the principal mo- 
ment of inertia. A nonaxisymmetric neutron star rotat- 
ing with frequency v emits periodic gravitational waves 
with amplitude 

/.o = ^^e, (2) 

where G is Newton's gravitational constant, c is the speed 
of light, /zz is the principal moment of inertia of the ob- 
ject, / (equal to ^v) is the gravitational wave frequency. 



and d is the distance to the object. Equation (jS)) gives 
the strain amplitude of a gravitational wave from an op- 
timally oriented source [see Eq. H25|) below]. 

The ellipticity of neutron stars is highly uncertain. The 
maximum ellipticity that can be supported by a neutron 
star's crust is estimated to be 

e„..«5xl0-^(^) , (3) 

where a is the breaking strain of the solid crust. The 
numerical coefficient in Eq. JSJ is small mainly because 
the shear modulus of the inner crust (which constitutes 
most of the crust's mass) is small, in the sense that it 
is about 10"'^ times the pressure. Eq. (O uses a fiducial 
breaking strain of 10~^ since that is roughly the upper 
limit for the best terrestrial alloys. However, a could be 
as high as 10""'^ for a perfect crystal with no defects [28| . 
or several orders of magnitude smaller for an amorphous 
solid or a crystal with many defects. 

Some exotic alternatives to standard neutron stars 
feature solid cores, which could support considerably 
larger cUipticitics . The most speculative and highcst- 
ellipticity model is that of a solid strange-quark star, for 
which 

.„.,»4xl0-'(^). (4) 

This much higher value of e^^^ is mostly due to the higher 
shear modulus, which for some strange star models can be 
almost as large as the pressure. Another (still speculative 
but more robust) model is the hybrid star, which consists 
of a normal neutron star outside a solid core of mixed 
quark and baryon matter, which may extend from the 
center to nearly the bottom of the crust. For hybrid 
stars, 

e^^gxlO-*^ (^), (5) 

although this is highly dependent on the poorly known 
range of densities occupied by the quark-baryon mixture. 
Stars with charged meson condensates could also have 
solid cores with overall cUipticitics similar to those of 
hybrid stars. 

Regardless of the maximum ellipticity supportable by 
shear stresses, there is the separate problem of how to 
reach the maximum. The crust of a young neutron star 
probably cracks as the neutron star spins down, but it is 
unclear how long it takes for gravity to smooth out the 
neutron star's shape. Accreting neutron stars in bina- 
ries have a natural way of reaching and maintaining the 
maximum deformation, since the accretion flow, guided 
by the neutron star's magnetic field, naturally produces 
"hot spots" on the surface, which can imprint themselves 
as lateral temperature variations throughout the crust. 
Through the temperature dependence of electron cap- 
ture, these variations can lead to "hills" in hotter areas 
which extend down to the dense inner crust, and with 
a reasonable temperature variation the ellipticity might 
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reach the maximum elastic value 0. The accreted ma- 
terial can also be held up in mountains on the surface by 
the magnetic field itself: The matter is a good conduc- 
tor, and thus it crosses field lines relatively slowly and 
can pile up in mountains larger than those supportable 
by elasticity alone Q . Depending on the field configura- 
tion, accretion rate, and temperature, the ellipticity from 
this mechanism could be up to lO""' even for ordinary 
neutron stars. 

Strong internal magnetic fields are another possible 
cause of ellipticity Differential rotation immediately 
after the core collapse in which a neutron star is formed 
can lead to an internal magnetic field with a large toroidal 
part. Dissipation tends to drive the symmetry axis of 
a toroidal field toward the star's equator, which is the 
orientation that maximizes the ellipticity. The resulting 
ellipticity is 



1.6 X 10-6 



• B ^ 
, 101= G/ 
B 



B < IQi^G, 



1.6 X 10-6 B>ioi5G, 



(6) 



where B is the root-mcan-square value of the toroidal 
part of the field averaged over the interior of the star. 
Note that this mechanism requires that the external field 
be much smaller than the internal field, since such strong 
external fields will spin a star out of the LIGO frequency 
band on a very short timescale. 

An alternative way of generating asymmetry is the r- 
modcs, fiuid oscillations dominated by the Coriolis restor- 
ing force. These modes may be unstable to growth 
through gravitational radiation reaction (the CFS insta- 
bility) under astrophysically realistic conditions. Rather 
than go into the many details of the physics and astro- 
physics, we refer the reader to a recent review of the 
literature and summarize here only what is directly rel- 
evant to our search: The r-modes have been proposed 
as a source of gravitational waves from newborn neutron 
stars and from rapidly accreting neutron stars P, 0| . 
The CFS instability of the r-modes in newborn neutron 
stars is probably not a good candidate for detection be- 
cause the emission is very short-lived, low amplitude, or 
both. Accreting neutron stars (or quark stars) are a bet- 
ter prospect for a detection of r-mode gravitational radi- 
ation because the emission may be long-lived with a duty 
cycle near unity [s^, |^ . 

Finally we consider free precession, i.e. the wobble of 
a neutron star whose symmetry axis does not coincide 
with its rotation axis. A large-amplitude wobble would 
produce 



hn 



10 
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1 kpc 



500 Hz 



(7) 



where 0^ is the wobble amplitude in radians. Such wob- 
ble may be longer lived than previously thought j^, but 
the amplitude is still small enough that such radiation 
is a target for second generation interferometers such as 
Advanced LIGO. 

In light of our current understanding of emission mech- 
anisms, the most likely sources of detectable gravitational 



waves are isolated neutron stars (through deformations) 
and accreting neutron stars in binaries (through defor- 
mations or r-modes). 



B. Isolated neutron stars 

The target population of this search is isolated rotat- 
ing compact stars that have not been observed electro- 
magnctically. Current models of stellar evolution sug- 
gest that our Galaxy contains of order 10^ neutron stars, 
while only of order 10^ are active pulsars. Up to now 
only about 1500 have been observed 113; there are nu- 
merous reasons for this, including selection effects and 
the fact that many have faint emission. Therefore the 
target population is a large fraction of the neutron stars 
in the Galaxy. 



1. Maximum expected signal amplitude at the Earth 

Despite this large target population and the variety 
of GW emission mechanisms that have been considered, 
one can make a robust argument, based on energet- 
ics and statistics, that the amplitude of the strongest 
gravitational-wave pulsar that one could reasonably hope 
to detect on Earth is bounded by /iq < 4 x 10"^'*. The ar- 
gument is a modification of an observation due to Bland- 
ford (which was unpublished, but credited to him in 
Thome's review in j33j). 

The argument begins by assuming, very optimistically, 
that all neutron stars in the Galaxy are born at very 
high spin rate and then spin down principally due to 
gravitational wave emission. For simplicity we shall also 
assume that all neutron stars follow the same spin-down 
law v{i') or equivalently /(/), although this turns out to 
be unnecessary to the conclusion. It is helpful to express 
the spin-down law in terms of the spin-down timescale 



.(/) 



/ 



|4/(/)| 



(8) 



For a neutron star with constant ellipticity, Tg„(/) is the 
time for the gravitational wave frequency to drift down 
to / from some initial, much higher spin frequency — 
but the argument does not place any requirements on 
the ellipticity or the emission mechanism. A source's 
gravitational wave amplitude is then related to Tgw (/) 

by 



gw 



(/)' 



(9) 



Here we are assuming that the star is not accreting, so 
that the angular momentum loss to GWs causes the star 
to slow down. The case of accreting neutron stars is dealt 
with separately, below. 
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We now consider the distribution of neutron stars 
in space and frequency. Let N{f)Af be the num- 
ber of Galactic neutron stars in the frequency range 
[/ - A//2, / + A//2]. We assume that the birthrate has 
been roughly constant over about the last 10^ years, so 
that this distribution has settled into a statistical steady 
state: dN{f)/dt = 0. Then N{f)f is just the neutron 
star birthrate 1/rb, where Tb may be as short as 30 years. 
For simplicity, we model the spatial distribution of neu- 
tron stars in our Galaxy as that of a uniform cylindrical 
disk, with radius Rq ~ lOkpc and height H « 600 pc. 
Then the density n(/) of neutron stars near the Earth, 
in the frequency range [/ — A//2,/ -I- A//2], is just 
n(/)A/ = [T,RlH)-^N{f)Af. 

Let N{f,r) be that portion of N{f) due to neutron 
stars whose distance from Earth is less than r. For H/2 < 
r < Rg, we have 



dNif,r) 
dr 



= 2TTrHn{f) 



(10) 
(11) 



(and it drops off rapidly for r > Rg)- Changing variables 
from r to using Eqs. (jSJ and we have 



dN{f,ho) 
dhn 



3 5GL 



2 c^TbRl 



(12) 



Note that the dependence on the poorly known Tgw{f) 
has dropped out of this equation. This was the essence 
of Blandford's observation. 

Now consider a search for GW pulsars in the frequency 
range [/„;„, /„,ax]- Integrating the distribution in Eq. H12fl 
over this band, we obtain the distribution of sources as a 
function of h^: 



dN, 



band 



dho 



h--^ In 



(13) 



The amplitude /iq of the strongest source is implicitly 
given by 



dNhi^nd 1 

^ 2 ■ 



(14) 



That is, even given our optimistic assumptions about 
the neutron star population, there is only a fifty percent 
chance of seeing a source as strong as /ig""'. The integral 
in Eq. (fTHl is trivial; it yields 



5G/z 



I 1/2 



In 



1/2 



(15) 



Inserting [ln(/„^^//„,i„)] « 1 (appropriate for a typical 
broadband search, as conducted here), and adopting as 
fiducial values /zz = lO^^gcm^, Rq = lOkpc, and t\, 
30 yr, we arrive at 



4 X 10 



-24 



(16) 



This is what we aimed to show. 

We now address the robustness of some assumptions 
in the argument. First, the assumption of a universal 
spin-down function Tg„(/) was unnecessary, since Tg„(/) 
disappeared from Eq. 1)12(1 and the subsequent equations 



that led to hJk 



Had we divided neutron stars into 



different classes labelled by i and assigned each a spin- 
down law Tg^(/) and birthrate l/r^, each would have 

contributed its own term to dN j dh^ which would have 
been independent of t"^^^ and the result for /iq'"" would 
have been the same. 

Second, in using Eq. (|10|1 . we have in effect assumed 
that the strongest source is in the distance range i//2 < 
r < Rq. We cannot evade the upper limit by assuming 
that the neutron stars have extremely long spin-down 
times (so that r < H/2) or extremely short ones (so that 
the brightest is outside our Galaxy, r > Rq). If the 
brightest sources are at r < H/2 (as happens if these 
sources have long spin-down times, Tgw > Tb{2RQ/H)^), 
then our estimate of /ig only decreases, because at short 
distances the spatial distribution of neutron stars be- 
comes approximately spherically symmetric instead of 
planar and the right hand sides of Eqs. (|10|l and l|12l) 
are multiplied by a factor 2r/H < 1. On the other 
hand, if Tgw(/) (in the LIGO range) is much shorter 
than Tb, then the probability that such an object ex- 
ists inside our Galaxy is ^ 1. For example, a neutron 
star with Tgw(/) = 3 yr located at r = lOkpc would 
have ho = 4.14 x 10~^^, but the probability of currently 
having a neutron star with this (or shorter) Tg„ is only 

^gw/Tb < 1/10. 

Third, we have implicitly assumed that each neutron 
star spins down only once. In fact, it is clear that some 
stars in binaries are "recycled" to higher spins by accre- 
tion, and then spin down again. This effectively increases 
the neutron star birth rate (since for our purposes the 
recycled stars are born twice), but since the fraction of 
stars recycled is very small the increase in the effective 
birth rate is also small. 



2. Expected sensitivity of the S2 search 

Typical noise levels of LIGO during the S2 run were 
approximately [Sh{f)Y^'^ « 3 x 10~'^^ Hz~^/^ where Sh 
is the strain noise power spectral density, as shown in 
Fig. ^ Even for a known GW pulsar with an average sky 
position, inclination angle, polarization, and frequency, 
the amplitude of the signal that we could detect in Gaus- 
sian stationary noise with a false alarm rate of 1% and a 
false dismissal rate of 10% is Il5l 



{hoif)) 



11.4a 



^obs 



(17) 



where Tobs is the integration time and the angled brackets 
indicate an average source. In all-sky searches for pulsars 
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with unknown parameters, the aniphtudc ho must be sev- 
eral times greater than this to rise convincingly above the 
background. Therefore, in Tabs = 10 hours of S2 data, 
signals with amplitude ho below about 10~^^ would not 
be detectable. This is a factor w 25 greater than the ft-jf 
of Eq. so our S2 analysis is unlikely to be sensitive 
enough to reveal previously unknown pulsars. 

The sensitivity of our search is further restricted by 
the template bank, which does not include the effects 
of signal spin-down for reasons of computational cost. 
Phase mismatch between the signal and matched filter 
causes the detection statistic (see Sec. IV A|) to decrease 
rapidly for GW frequency derivatives / that exceed 

-ax[/] ^ 1t-J = 4 X 10-1° 'hzs"1 . (18) 

Assuming that all of the spin-down of a neutron star 
is due to gravitational waves (from a mass quadrupole 
deformation), our search is restricted to pulsars with el- 
lipticity e less than 

/ 5c^max[/] 

This limit, derived from combining the quadrupole for- 
mula for GW luminosity 

(the first factor is 1/10 instead of 1/5 due to time aver- 
aging of the signal) with the kinetic energy of rotation 

i? = i^V'4z, (21) 
(assuming / = 2i^) takes the numerical value 

.,.c,^eixlO-(^— j (_) (22) 

for our maximum /. 

The curves in Fig.|21are obtained by combining Eqs. Q 
and H17|l "^ and solving for the distance d for different val- 
ues of the ellipticity, using an average value for noise in 
the detectors during the S2 run. The curves show the 



^ Note that the value of hg derived from Eq. 1171 yields a value 
of the detection statistic 2jF for an average source as seen with 
a detector at S2 sensitivity and over an observation time of 10 
hours, of about 21, which is extremely close to the value of 20 
which is used in this analysis as threshold for registering candi- 
date events. Thus combining Eqs. J3J and IIYI determines the 
smallest amplitude that our search pipeline could detect (corre- 
sponding to a signal just at the threshold), provided appropriate 
follow-up studies of the registered events ensued. 



70 




Frequency (Hz) 

FIG. 2: Effective average range (defined in the text) of our 
search as a function of frequency for three ellipticities: 10~^ 
(maximum for a normal neutron star), 10^'' (maximum for a 
more optimistic object), and esd, the spin-down limit defined 
in the text. Note that for sources above 300 Hz the reach of 
the search is limited by the maximum spin-down value of a 
signal that may be detected without loss of sensitivity. 

average distance, in the sense of the definition H17|l . at 
which a source may be detected. 

The dark gray region shows that a GW pulsar with 
e = 10^^ could be detected by this search only if it 
were very close, less than ~ 5 parsecs away. The light 
gray region shows the distance at which a GW pulsar 
with e = 10~^ could be detected if templates with suffi- 
ciently large spin-down values were searched. However, 
this search can detect such pulsars only below 300 Hz, 
because above 300 Hz a GW pulsar with e = 10"^ spins 
down too fast to be detected with the no-spin-down tem- 
plates used. The thick line indicates the distance limit 
for the (frequency-dependent) maximum value of epsilon 
that could be detected with the templates used in this 
search. At certain frequencies below 300 Hz, a GW pul- 
sar could be seen somewhat farther away than 30 pc, but 
only if it has e > 10~^. Although egd and the corre- 
sponding curve were derived assuming a quadrupolar de- 
formation as the emission mechanism, the results would 
be similar for other mechanisms. Equation H21|l includes 
an implicit factor /^/(2i^)^, which results in esd and the 
corresponding range (for a fixed GW frequency /) being 
multiplied by f /{2i/), which is 1/2 for free precession and 
about 2/3 for 7'-modes. Even for a source with optimum 
inclination angle and polarization, the range increases 
only by a factor « 2. The distance to the closest known 
pulsar in the LIGO frequency band, PSR J0437-4715, 
is about 140 pc The distance to the closest known 
neutron star, RX J185635-3754, is about 120 pc [H. 
Therefore this search would be sensitive only to nearby 
previously unknown objects. 

While we have argued that a detection would be very 
unlikely, it should be recalled that Eq. Hlfci|) was based on 
a statistical argument. It is always possible that there is 
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a GW-bright neutron star that is much closer to us than 
would be expected from a random distribution of super- 
novae (for example due to recent star formation in the 
Gould belt as considered in Is^ ) . It is also possible that a 
"blind" search of the sort performed here could discover 
some previously unknown class of compact objects not 
born in supernovae. 

More importantly, future searches for previously undis- 
covered rotating neutron stars using the methods pre- 
sented here will be much more sensitive. The goal of 
initial LIGO is to take a year of data at design sensitiv- 
ity. With respect to S2, this is a factor 10 improvement 
in the amplitude strain noise at most frequencies. The 
greater length of the data set will also increase the sen- 
sitivity to pulsars by a factor of a few (the precise value 
depends on the combination of coherent and incoherent 
analysis methods used). The net result is that initial 
LIGO will have Hq reduced from the S2 value by a factor 
of 30 or more to a value comparable to /ig 4 x 10^^"' 
of Eq. (Cni). 

C. Accreting neutron stars 



cal maximum spin frequencies, estimated at ~ 1400 Hz, 
the observed spin distribution is hard to explain without 
some competing mechanism, such as gravitational radia- 
tion, to halt the spin-up. Since the gravitational torque 
scales as i^^, gravitational radiation is also a natural ex- 
planation for why the spin frequencies occupy a rather 
narrow window: a factor 32 difference in accretion rate 
leads to only a factor 2 difference in equilibrium spin 
rate Q. 

If the above argument holds, then the accreting neu- 
tron star brightest in X-rays is also the brightest in grav- 
itational waves. Sco X-1, which was the first extrasolar 
X-ray source discovered, is the strongest persistent X-ray 
source in the sky. Assuming equilibrium between GWs 
and accretion, the gravitational wave strain of Sco X-1 
at the Earth is 

,..3xl0-»(^)"\ (24, 

which should be detectable by second generation inter- 
ferometers. The gravitational wave strains from other 
accreting neutron stars are expected to be lower. 



Maximum expected signal amplitude at Earth 



2. Expected sensitivity of S2 search for Sco X-1 



The robust upper limit in Eq. H16() refers only to non- 
accreting neutron stars, since energy conservation plays 
a crucial role. If accretion replenishes the star's angu- 
lar momentum, a different but equally robust argument 
(i.e., practically independent of the details of the emis- 
sion mechanism) can be made regarding the maximum 
strain /ig'"^ at the Earth. In this case /iq is set by the 
X-ray luminosity of the brightest X-ray source. 

The basic idea is that if the energy (or angular momen- 
tum) lost to GWs is replenished by accretion, then the 
strongest GW emitters are those accreting at the highest 
rate, near the Eddington limit. Such systems exist: the 
low-mass X-ray binaries (LMXBs) , so-called since the ac- 
creted material is tidally stripped from a low-mass com- 
panion star. The accreted gas hitting the surface of the 
neutron star is heated to 10* K and emits X-rays. As 
noted several times over the years 0, [s^ , if one as- 
sumes that spin-down from GW emission is in equilib- 
rium with accretion torque, then the GW amplitude /ig 
is directly related to the X-ray luminosity: 



5 X 10 



-27 



300 Hz 



1/2 / p \ 1/2 



10 = erg cm 



(23) 

where is the X-ray flux. In the 1970s when this con- 
nection was first proposed, there was no observational 
support for the idea that the LMXBs are strong GW 
emitters. But the spin frequencies of many LMXBs 
are now known, and most are observed to cluster in a 
fairly narrow range of spin frequencies 270 Hz ^ ^ 
620 Hz [s^. Since most neutron stars will have accreted 
enough matter to spin them up to near their theoreti- 



The orbital parameters of Sco X-1 are poorly con- 
strained by present (mainly optical) observations and 
large uncertainties affect the determination of the rota- 
tion frequency of the source (details are provided in Sec- 
tion ^^^^l- The immediate implication for a coherent 
search for gravitational waves from such a neutron star 
is that a very large number of discrete templates are re- 
quired to cover the relevant parameter space, which in 
turn dramatically increases the computational costs |39l| . 
The optimal sensitivity that can be achieved with a co- 
herent search is therefore set primarily by the length of 
the data set that one can afford to process (with fixed 
computational resources) and the spectral density of the 
detector noise. As we discuss in Section HV B 21 the max- 
imum span of the observation time set by the computa- 
tional burden of the Sco X-1 pipeline (approximately one 
week on « 100 CPUs) limits the observation span to 6 
hours. 

The overall sensitivity of the search that we are de- 
scribing is determined by each stage of the pipeline, 
which we describe in detail in Section IV Bl Assuming 
that the noise in the instrument can be described as a 
Gaussian and stationary process (an assumption which 
however breaks down in some frequency regions and/or 
for portions of the observation time) we can statistically 
model the effects of each step of the analysis and es- 
timate the sensitivity of the search. The results of such 
modelling through the use of Monte Carlo simulations are 
shown in Fig |21 where we give the expected upper limit 
sensitivity of the search implemented for the analysis. We 
contrast this with the hypothetical case in which the Sco 
X-1 parameters are known perfectly making it a single 



9 



1 

Jl 1 I ill 


II 

k 1 1 


1 

^^jjLJ>Aj- — ■ — J—iL 











200 300 400 500 600 700 800 900 1000 

Frequency (Hz) 



FIG. 3: Here we show the expected upper hmit sensitivity of 
the S2 Sco X-1 search. The upper black curve represents the 
expected sensitivity of the S2 analysis based on an optimally 
selected 6 hr dataset (chosen specifically for our search band). 
The gray curve (second from the top) shows the sensitivity 
in the hypothetical case in which all of the Sco X-1 system 
parameters are known exactly making Sco X-1 a single filter 
target and the entire S2 data set is analyzed. Both curves are 
based on a 95% confidence upper limit. The remaining curves 
represent V Sh (/) /Tobs for LI (black) and HI (gray); Sh{f) 
is the typical noise spectral density that characterizes the LI 
and HI data, and Tobs is the actual observation time (taking 
into account the duty cycle, which is different for LI and HI) 
for each instrument. 



filter target for the whole duration of the S2 run. The 
dramatic difference (of at least an order of magnitude) 
between the estimated sensitivity curves of these two sce- 
narios is primarily due to the large parameter space we 
have to search. This has two consequences, which con- 
tribute to degrading the sensitivity of the analysis: (i) we 
are computationally limited by the vast number of tem- 
plates that we must search and therefore must reduce the 
observation to a subsection of the S2 data, and (ii) sam- 
pling a large number of independent locations increases 
the probability that noise alone will produce a high value 
of the detection statistic. 



We note that the S2 Sco X-1 analysis (see Sec- 
tion ^^^^ is a factor of « 5000 less sensitive than the 
characteristic amplitude given in Eq. (HH). In the hy- 
pothetical case in which Sco X-1 is a single filter target 
and we are able to analyze the entirety of S2 data, then 
we are still a factor ^ 100 away. However, as mentioned 
in the introduction, the search reported in this paper 
will be one of the stages of a more sensitive "hierarchi- 
cal pipeline" that will allow us to achieve quasi-optimal 
sensitivity with fixed computational resources. 



IV. SIGNAL MODEL 
A. The signal at the detector 

We consider a rotating neutron star with equatorial co- 
ordinates a (right ascension) and 6 (declination). Gravi- 
tational waves propagate in the direction k and the star 
spins around an axis whose direction, assumed to be con- 
stant, is identified by the unit vector s. 

The strain h(t) recorded at the interferometer output 
at detector time t is: 

h{t) ^ ho ^ (l + cos^t)F+(i;a,(5,i/') cos$(t) 

+ cos LFx{t] a, 5, ip) sin^{t) , (25) 

where V is the polarization angle, defined as tantjj = 
[{s ■ k) (z ■ k) — {s ■ k)]/k ■ {s X z), z is the direction to 
the north celestial pole, and cost ~ k ■ s. Gravitational 
wave laser interferometers are all-sky monitors with a 
response that depends on the source location in the sky 
and the wave polarization: this is encoded in the (time 
dependent) antenna beam patterns -F+^x (i; ck, The 
term $(i) in Eq. (|25|) represents the phase of the received 
gravitational signal. 

The analysis challenge to detect weak quasi-periodic 
continuous gravitational waves stems from the Doppler 
shift of the gravitational phase ^{t) due to the relative 
motion between the detector and the source. It is con- 
venient to introduce the following times: t, the time 
measured at the detector; T, the solar-system-barycentcr 
(SSB) coordinate time; and tp, the proper time in the rest 
frame of the pulsar^. 

The timing model that links the detector time t to the 
coordinate time T at the SSB is: 



T = t + 



iS0 : 



(26) 



where r is a (time-dependent) vector from the SSB to 
the detector at the time of the observations, ti is a unit 
vector towards the pulsar (it identifies the source posi- 
tion in the sky) and A^© and Ag© are the solar system 
Einstein and Shapiro time delays, respectively |43|. For 
an isolated neutron star tp and T arc equivalent up to 
an additive constant. If the source is in a binary sys- 
tem, as it is the case for Sco X-1, significant additional 
accelerations are involved, and a further transformation 
is required to relate the proper time tp to the detector 
time t. Following ji^l, we have: 



T-Tn = U 



Af 



A, 



(27) 



^ Notice that our notation for the three different times is differ- 
ent from the est ablis hed conventions adopted in the radio pulsar 
community, e.g. l4Cl . 
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where A/j, the Roemer time delay, is analogous to the 
solar system term (r ■ n) /c; A^; and A5 arc the orbital 
Einstein and Shapiro time delay, analogous to Aeq and 
Asq; and Tq is an arbitrary (constant) reference epoch. 
For the case of Sco X-1, we consider a circular orbit 
for the analysis (cf Section FlV B 21 for more details) and 
therefore set Ae = 0- Furthermore, the binary is non- 
relativistic and from the source parameters we estimate 
A5 < 3 /iS which is negligible. For a circular orbit, the 
Rocmer time delay is simply given by 

Ar = ^ sin(u + uj) (28) 
c 

where Op is the radius of the neutron star orbit projected 
on the line of sight, u the argument of the periapsis and u 
the so-called eccentric anomaly; for the case of a circular 
orbit u = 2'!T{tp — tp,Q) / P, where P is the period of the bi- 
nary and tp^o is a constant reference time, conventionally 
referred to as the "time of periapse passage" . 

In this paper we consider gravitational waves whose 
intrinsic frequency drift is negligible over the integration 
time of the searches (details are provided in the next 
section), both for the blind analysis of unknown isolated 
neutron stars and Sco X-1. The phase model is simplest 
in this case and given by: 

$(ip) = 2^/otj, + $o, (29) 

where $0 is an overall constant phase term and /o is the 
frequency of the gravitational wave at the reference time. 

B. Parameter space of the search 

The analysis approach presented in this paper is used 
for two different searches with different search parame- 
ters. Both searches require exploring a three dimensional 
parameter space, made up of two "position parameters" 
(whose nature is different for the two searches) and the 
unknown frequency of the signal. For the all-sky blind 
analysis aimed at unknown isolated neutron stars one 
needs to Dopplcr correct the phase of the signal for any 
given point in the sky, based on the angular resolution of 
the instrument over the observation time, and so a search 
is performed on the sky coordinates a and i5. For the Sco 
X-1 analysis, the sky location of the system is known, 
however the system is in a binary orbit with poorly mea- 
sured orbital elements; thus, one needs to search over a 
range of orbital parameter values. The frequency search 
parameter is for both searches the /o defined by Eq. (|29|l . 
where the reference time has been chosen to be the time- 
stamp of the first sample of the data set. The frequency 
band over which the two analyses are carried out is also 
different, and the choice is determined by astrophysical 
and practical reasons. As explained in Sec. IV CI the data 
set in HI does not coincide in time with the LI data set 
for either of the analyses. Consequently a signal with a 
non-zero frequency derivative would appear at a differ- 
ent frequency template in each data set. However, for 



the maximum spin-down rates considered in this search, 
and given the time lag between the two data sets, the 
maximum difference between the search frequencies hap- 
pens for the isolated objects search and amounts to 0.5 
mHz. We will see that the frequency coincidence win- 
dow is much larger than this and that when we discuss 
spectral features in the noise of the data and locate them 
based on template-triggers at a frequency /o, the spec- 
tral resolution is never finer than 0.5 mHz. So for the 
practical purposes of the present discussion we can ne- 
glect this difference and will often refer to /o generically 
as the signal's frequency. 



1. Isolated neutron stars 

The analysis for isolated neutron stars covers the entire 
sky and we have restricted the search to the frequency 
range 160-728.8 Hz. The low frequency end of the band 
was chosen because the depth of our search degrades sig- 
nificantly below 160 Hz, sec Fig. |2 The choice of the 
high frequency limit at 728.8 is primarily determined by 
the computational burden of the analysis, which scales 
as the square of the maximum frequency that is searched 
for. 

In order to keep the computational costs at a reason- 
able level (< 1 month on ^ 800 CPUs), no explicit 
search over spin-down parameters was carried out. The 
length of the data set that is analyzed is approximately 
10 hours, thus no loss of sensitivity is incurred for sources 
with spin-down rates smaller than 4 x 10^^" Hzs~^; see 
Eq. H18|) . This is a fairly high spin-down rate compared 
to those measured in isolated radio pulsars; however it 
does constrain the sensitivity for sources above 300 Hz, 
as can be seen from Fig. |21 



2. Sco X-1 

Sco X-1 is a neutron star in a 18.9 h orbit around a low 
mass (~ 0.42 Mq) companion at a distance r — 2.8 ± 0.3 
kpc from Earth. In this section we review our present 
knowledge of the source parameters that are relevant for 
gravitational wave observations. Table ^contains a sum- 
mary of the parameters and the associated uncertainties 
that define the search area. In what follows wc will as- 
sume the observation time to be 6 hours. This is approx- 
imately what was adopted for the analysis presented in 
this paper. We will justify this choice at the end of the 
section. 

The most accurate determination of the Sco X-1 sky 
position comes from Very Long Baseline Array (VLBA) 
observations 0,E3 is reported in Tabled The over- 
all error on the source location is ~ 0.5 arcsec, which is 
significantly smaller than the 100 arcsec sky resolution 
associated with a 6 hour GW search. Hence we assume 
the position of Sco X-1 (i.e. the barycenter of the binary 
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right ascension 




16h 19m 55.0850s 


declination 


s 


— 15° 38' 24 9" 


proper motion (east-west direction) 




—0.00688 ± 0.00007 arcsec yr~ 


proper motion (north-south direction) 




0.01202 ± 0.00016 arcsec yr~^ 


distance 


d 


2.8 ± 0.3 kpc 


orbital period 


P 


68023.84 ± 0.08 sec 


time of periapse passage 


T 


731163327 ± 299 sec 


projected semi-major axis 


Up 


1.44 ±0.18 sec 


eccentricity 


e 


< 3 X 10"^ 


QPOs frequency separation 


237 ± 5 Hz < Az^QPo < 307 ± 5 H: 



TABLE I: The parameters of the low-mass X-ray binary Scorpius X-1. The quoted measurement errors are all l-cr. We refer 
the reader to the text for details and references. 



system) to be exactly known and we "point" (in software) 
at that region of the sky. 

Three parameters describe the circular orbit of a star 
in a binary system: the orbital period (P), the projection 
of the semi-major axis of the orbit Op (which for e = 
corresponds to the projected radius of the orbit) and the 
location of the star on the orbit at some given reference 
time. For eccentric orbits this is usually parameterized 
by the time of periapse passage (or time of periastron). 
In the case of a circular orbit we define the orbital phase 
reference time T as the time at which the star crosses the 
ascending node as measured by an observer at the SSB. 
This is equivalent to setting the argument of periapse 
(the angle between the ascending node and the direction 
of periapsis) to zero. 

In the case of Sco X-1, P is by far the most accurately 
determined parameter [43j, and over a 6 hour search it 
can be considered known because the loss of signal-to- 
noise ratio (SNR) introduced by matching two templates 
with any value of P in the range of Table is negligible. 
P becomes a search parameter, requiring multiple filters, 
only for coherent integration times ^ 10^ s. The major 
orbital parameters with the largest uncertainties are the 
projected semi-major axis of the orbit along the line of 
sight, Op, and the orbital phase reference time. The large 
uncertainty on Cp is primarily due to the poor determi- 
nation of the orbital velocity (40 ± 5kms~^ Q). The 
uncertainty on the orbital phase reference time is due to 
the difficulty in locating the Sco X-1 low-mass compan- 
ion on the orbit. The search therefore requires a discrete 
grid of filters in the (opjT) space. 

We assume that Sco X-1 is in a circular orbit, which 
is what one expects for a semi-detached binary system 
and which is consistent with the best fits of the orbital 
parameters However, orbital fits for models with 

e ^ were clearly dominated by the noise introduced by 
the geometry of the Roche lobe Over an integra- 

tion time of ~ 6 h, the eccentricity needs to be smaller 
than ~ 10^"' in order for the detection statistic !F to be 
affected less than 1%; for e ~ 10~^, losses of the order 
of 10% are expected and are consistent with the results 



presented later in the paper. Unfortunately current ob- 
servations are not able to constrain e to such levels of 
accuracy: in this paper we adopt the strategy of analyz- 
ing the data under the assumption e = 0, and we quantify 
(for a smaller set of the parameter space) the efficiency of 
the pipeline in searching for gravitational waves emitted 
by a binary with non zero eccentricity; in other words, we 
quote upper limits for different values of the eccentricity 
that are obtained with non-optimal search templates. 

The last parameter we need to search for is the fre- 
quency of the gravitational radiation /. The rotation 
frequency v of Sco X-1 is inferred from the difference 
of the frequency of the kHz quasi periodic oscillations 
(QPOs). Unfortunately this frequency difference is not 
constant, and over a 4 day observation has shown a 
very pronounced drift between 237 ± 5 Hz to 307 ± 5 Hz, 
where the errors should be interpreted as the la val- 
ues This drifting of QPO frequency separation was 
found to be positively correlated to the inferred mass ac- 
cretion rate. It is also important to stress that there is 
a still unresolved controversy as to whether the adopted 
model that links v to the difference of the frequency of 
the kHz QPOs is indeed the correct one, and if it is valid 
for all the observed LMXBs. Moreover, the gravitational 
wave frequency / is related to in a different way, de- 
pending on the model that is considered: f = 2v \i one 
considers nonaxisymmetric distortions and / = (4/3)z/ 
if one considers unstable r-modes. It is therefore clear 
that a search for gravitational waves from Sco X-1 should 
assume that the frequency is essentially unknown and 
the whole LIGO sensitivity band (say from 100 Hz to 
K, 1 kHz) should be considered. Because of the heavy 
computational burden, such a search requires a differ- 
ent approach (this search is currently in progress). For 
the analysis presented in this paper, we have decided to 
confine the search to GWs emitted by nonaxisymmet- 
ric distortions (/ = 2;^) and to constrain the frequency 
band to the two 20 Hz wide bands (464-484 Hz and 604- 
624 Hz) that bound the range of the drift of i/, according 
to currently accepted models for the kHz QPOs. The to- 
tal computational time for the analysis can be split into 
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two parts: (i) the search time Tscarch needed to search the 
data and, if no signal is detected, (ii) the upper limit time 
Tinj required to repeatedly inject and search for artifi- 
cially generated signals for the purposes of setting the up- 
per limits. Let Tgpan be the span of the data set which is 
analyzed, that is the difference between the time stamps 
of the first and last data point in the time series. Let Tobs 
be the effective duration of the data set containing non- 
zero data points. The definitions imply Tobs < ?span, and 
for data with no gaps Tobs = ?span- For a search confined 
to a period (sufficiently) shorter than the orbital period 
of the source, the two computational times are: 




where A/ is the search frequency band, AT the search 
range for the time of periapse passage, fj, is the template 
bank mismatch, and A'cpu is the number of ~ 2 GHz 
CPUs available The quantities iVtiiais and Nh^ are 
the number of artificial signals injected per value of Hq 
and the number of different values of ho injected, respec- 
tively. Note the steep dependency of the search time 
Tscaich on the maximum observation time span Tspan- 
The contributing factors to this scaling are the increas- 
ing number of orbital and frequency filters, iVorb and 
-^froq respectively, with observation time span, where 
iVorb oc Tgpan and A'froq OC Tobs. Thcrc is also a linear 
scaling of computational time with Tgpan (corrected by 
the factor Tobs /Tspan that takes into account only the 
non-zero data points) due to increased data volume be- 
ing analyzed. Prom Eqs. I|3U|I and H31|) it is therefore 
clear that if one wants to complete the full analysis over 
a period ^ 1 week the choice Tspan = 6 h is appropriate. 

V. ANALYSIS OF THE DATA 

In this section we describe the analysis strategy and its 
implementation on the data collected during the second 
science run by the two 4-km LIGO interferometers. 

The inner core of the analysis is built on the frequency- 
domain matched-filter approach that we applied to the 
data collected during the first science run to place 
an upper limit on gravitational radiation from PSR 
J1939-|~2134 However, this analysis is considerably 
more complex with respect to |l5j because (i) the search 
is carried out over a large number of templates (either 



over sky position or source orbital parameters), (ii) the 
data are analyzed in coincidence between two interferom- 
eters in order to reduce the false alarm probability and 
thereby improve the overall sensitivity of the search, and 
(iii) the upper limit is derived from the maximum joint 
significance of coincident templates. 

A. The detection statistic 

The optimal detection statistic (in the maximum 
likelihood sense) to search coherently for quasi- 
monochromatic signals is the so-called ^-statistic'^ in- 
troduced in This statistic can be extended in a 
straightforward manner to the case of a signal from a 
pulsar in a binary system. 

In the absence of signal, is distributed according 
to a (central) distribution with 4 degrees of freedom 
and the relevant probability density function is given by 

Po(2^) = . (32) 

We define the false alarm probability of as 

P^{2T)^ / po{2T')d{2T'). (33) 

In the presence of a signal, 2T follows a non-central 
distribution with 4 degrees of freedom and non-centrality 
parameter p^; the associated probability density function 
is 

p^{2T) = ie-(2^+''^)/^ & (v/277) , (34) 
2 VP 

where Ii is the modified Bessel function of the first kind 
of order one and 

^ SIT) f" <'^' 

The expected value of 2T \s A + p^. Prom Eq. it is 
clear that the detection statistic is proportional to the 
square of the amplitude of the gravitational wave signal, 
hi, given by Eq. (j2Jl. 

B. Pipeline 

The search pipeline is schematically illustrated in 
Pig.^ A template bank is set up for each search covering 
the parameter space under inspection. Por both analyses 



^ We would like to stress that this statistic is completely unrelated 
to the _F-statistic described in statistical textbooks to test the 
null hypothesis for two variances drawn from distributions with 
the same mean. 
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the template bank is three-dimensional: it covers right 
ascension and declination for the unknown isolated pul- 
sar search, and the orbital phase reference time and the 
projection of the orbital semi-major axis for the Sco X- 
1 analysis. In addition in both cases we search for the 
unknown gravitational wave frequency. 

The data stream is treated in exactly the same way for 
each search: the full search frequency band is divided into 
smaller 1 Hz) sub-bands'^, the J^-statistic is computed 
at every point in the template bank, and lists of candi- 
date templates are produced. Search template values are 
recorded when the detection statistic exceeds the value 
2J^ = 20, and we will refer to them as registered tem- 
plates. Note that we will also refer to these templates 
as "events", by analogy with the time-domain matched 
filtering analysis. 

In the frequency domain search conducted in [l5| we 
searched a single template in four detectors. Here we 
search a total of 5 x 10^^ templates in each detector 
for the isolated pulsars and 3 x 10^° templates overall 
in two detectors for Sco X-1. In order to reduce the 
number of recorded templates, only the maximum of the 
detection statistic over a frequency interval at fixed val- 
ues of the remaining template parameters is stored. The 
frequency interval over which this maximization is per- 
formed is based on the maximum expected width of the 
detection statistic for an actual signal. 

For each recorded template from one detector, the list 
of templates from the other detector is scanned for tem- 
plate(s) close enough so that an astrophysical signal could 
have been detected in both templates. The criteria used 
to define "closeness" are different for the two searches 
and will be described in sections IV F II and IVF 21 This 
procedure yields a third list of templates that are what 
we refer to as the coincident templates. These are tem- 
plate values for which the JT-statistic is above threshold 
and such that they could be ascribed to the same physi- 
cal signal in both data streams. The coincident templates 
are then ranked according to their joint significance. 

The most significant coincident template is identified in 
each ^ 1 Hz sub-band; we also refer to this as the loudest 
event for that frequency sub-band. An upper limit on 
the value of ho from a population of isolated sources, or 
from a family of possible source parameters in the Sco 
X-1 system, is placed in each frequency sub-band based 
on its loudest coincident event. Following this is 
done by injecting in the real data a set of fake signals 
at the same level of significance as the loudest measured 
event and by searching the data with the same pipeline 



Each of the sub-bands corresponds to the frequency region which 
the loudest candidate over the sky or the orbital parameters is 
maximized over and should be of comparable width with respect 
to typical noise floor variations. The precise bandwidth size is 
dictated by convenience in the computational set-up of the anal- 
yses, reflecting a different distribution of the computational load 
among the various nodes of the computer clusters used. 
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FIG. 4: Workflow of the pipeline. 



as was used in the analysis. The upper limit procedure 
is described in Sec. IV Gl 



C. Selection of the data set 

The data input to the search is in the form of short time 
baseline Fourier transforms of the time-domain data. At 
fixed observation time the computational cost of a search 
increases linearly with the number of short Fourier trans- 
forms (SFTs) employed. Hence, the longer the time base- 
line of the SFT, the less computationally intensive the 
search. There are two constraints to making the SFT 
time baseline long: i) the noise is estimated on the SFT 
time scale and thus it should be reasonably stationary on 
such timescale; and ii) the signal-to-noise ratio of a puta- 
tive source will be significantly degraded if the Doppler 
modulation during the SFT time baseline is of the order 
of 1/TsFT- For the S2 data set and a search extending to 
about 750 Hz we chose 30 minutes as the time baseline for 
the SFTs of the search for signals from isolated sources. 
The Sco X-1 search was carried out using 60 s long SFTs 
due to the more significant acceleration produced by the 
orbital motion. 

As described in the SFT data is normalized by the 
noise spectral amplitude. This quantity is estimated for 
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each SFT from the actual data near to the frequency bin 
of interest. In we used a simple average over frequen- 
cies around the target search frequency. That approach 
worked well because in the vicinity of the target fre- 
quency there were no spectral disturbances. But clearly 
we cannot count on this being the case while searching 
over several hundred Hz. We have thus adopted a spec- 
tral running median estimate method [SCI ,51^, 32] ■ In the 
isolated pulsar search we have chosen a very conservative 
window size of 50 1800 s-baseline-SFT bins (27.8 mHz) 
corresponding to a little under twice the number of terms 
used in the demodulation routine that computes the de- 
tection statistic through the integrals (108) and (109) of 
[4^ . We estimate the noise at every bin as the median 
computed on 25+25+1 values, corresponding to the 25 
preceding bins, the bin itself, and the 25 following bins. 
If an outlier in the data were due to a signal, our spectral 
estimate would be insensitive to it, and thus we would be 
preserving it in the normalized data. A window size of 
50 60 s-baseline-SFT frequency bins (0.833 Hz) was also 
used for the upper frequency band of the Sco X-1 search, 
604-624 Hz. Due to the presence of some large spectral 
features in the lower band, 464-484 Hz, a window size of 
25 60 s-baseline-SFT bins (0.417 Hz) was used in an at- 
tempt to better track the noise floor. Noise disturbances 
are evident in Fig. |3 where we show (with frequency 
resolution 1 /60 Hz) the average noise spectral density of 
the data set used in the analysis. Notice that the lower 
frequency band presents numerous spectral features, es- 
pecially in HI; moreover a strong and broad (approxi- 
mately 2 Hz) excess noise in both detectors is evident 
around 480 Hz, which corresponds to a harmonic of the 
60 Hz power line. 

The reconstruction of the strain from the output of the 
interferometer is referred to as the calibration. Details re- 
garding the calibration for the S2 run can be found in [23| • 
Both analyses presented here use a calibration performed 
in the frequency domain on SFTs of the detector output. 
The SFT-strain h{f) is computed by constructing a re- 
sponse function R(f,t) that acts on the interferometer 
output q{f): h{f) = R{f,t)q{f). Due mainly to changes 
in the amount of light stored in the Fabry-Perot cavi- 
ties of the interferometers, the response function, R{ f, t), 
varies in time. These variations arc measured using sinu- 
soidal excitations injected into the instrument. Through- 
out S2, changes in the response were computed every 60 
seconds. The SFTs used were 30 minutes long for the 
isolated pulsar analysis and an averaging procedure was 
used to estimate the response function for each SFT. For 
the binary search, which uses 60 s SFTs, a linear inter- 
polation was used, since the start times of the SFTs do 
not necessarily correspond to those at which the changes 
in the response were measured. 

The observation time chosen for the two searches is 
significantly less than the total observation time of S2, 
due to computational cost constraints: about 10 hours 
and 6 hours for the isolated pulsar and Sco X-1 searches, 
respectively. We picked the most sensitive data stretches 
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FIG. 5: The average of the noise over various 1 Hz sub-bands 
as described in the text for different sets of 20 SFTs from data 
of the LI detector (top plot) and HI detector (bottom plot). 
The X-axis labels the order number of the first SFT in each 
set. SFT #1 is the first SFT of the run. Neighboring sets only 
differ by one SFT. The cross indicates the data set chosen for 
the search for signals from isolated objects. 

covering the chosen observation times; the criteria used 
to select the data sets are described below and the dif- 
ferences reflect the different nature of the searches. 



1. Data selection for the isolated neutron star search 

Since the blind search for isolated neutron stars is an 
all-sky search, the most sensitive data are chosen based 
only on the noise performance of the detectors. The sen- 
sitivity is evaluated as an average of the sensitivity at 
different frequencies in the highest sensitivity band of 
the instrument. In particular, the noise is computed in 
six sub-bands that span the lowest 300 Hz range to be 
analyzed. The sub-bands are 1 Hz wide, with lowest fre- 
quencies respectively at 162 Hz, 219 Hz, 282 Hz, 338 Hz, 
398 Hz and 470 Hz. These sub-bands were chosen in re- 
gions free of spectral disturbances and the average power 
in these frequency regions can be taken as a measure of 
the noise floor there. Even though the search band ex- 
tends up to 730 Hz, we have chosen these reference sub- 
bands to lie below 500 Hz, because this is the most sen- 
sitive frequency range of our instruments. We construct 
sets of 20 SFTs (10 hrs of data), with the constraint that 
the data employed in each set does not span more than 
13 hours. This constraint stems from computational re- 
quirements: the spacing used for the template grid in the 
sky shrinks very fast with increasing spanned observation 
time. If the data contains no gaps then each 10 hr set 
differs from the previous only by a single SFT. For HI we 
arc able to construct 892 such sets, for LI only 8, reflect- 
ing the rather different duty cycle in the two instruments. 
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FIG. 6: The average amplitude spectral density \/Sh{f) of 
the data of the two detectors used for the isolated pulsars 
analysis. The bottom plot is the average over 1.2 Hz wide 
sub-bands of the average of the top two plots. The frequency 
resolution of the top two plots is Hz. The frequency 

resolution of the bottom plots is 2160 times coarser. 



This is obvious from the plots of Fig. |SJ For HI we were 
able to cover with sets of 20 SFTs the entire run in a 
fairly uniform way. For LI it was possible to find sets 
of nearly-contiguous SFTs only in the first and second 
quarter of the run. We finally compute the average over 
the different frequency sub-bands and wc pick the set for 
which this number is the smallest. Figure shows this 
average over the frequency bands and the cross points to 
the lowest- noise SFT-set. 

The data sets chosen were: for HI 20 30-minute SFTs 
starting at GPS time 733803157 that span 10 hours, 
and for LI 20 30-minute SFTs starting at GPS time 
732489168 spanning 12.75 hours. The plots of Fig. 
show the average power spectral density of this data set 
for the two detectors separately (top two plots) and the 
average of the same data over 1.2 Hz wide sub-bands and 
over the two detectors (bottom plot). 



2. Data selection for the Sco X-1 search 

We choose to analyze in each detector the most sensi- 
tive S2 data set which does not span more than 6 hours, 
which we have fixed based on computational cost con- 
straints. To rank the sensitivity of a data set we use the 
figure of merit 



Q{T) = 



HSh{T)) 



\ A{f) + B{f) 



^obs 



(36) 



where T = {Tstart, Tobs, Japan} identifies the data set — 
the time of the first data point Tgtart, the effective time 
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FIG. 7: The amplitude spectral density \/ Sh (/) of the opti- 
mally chosen data sets for LI and HI and for both frequency 
bands averaged over each 60 sec SFT. The solid black and 
dashed gray lines correspond to HI and LI, respectively. Note 
that the lower band contains a feature common to both de- 
tectors, a 60 Hz power line harmonic at 480 Hz with a width 
of ~ 2 Hz. The HI data set also contains a variety of other 
features in the lower band, some equally as large as the power 
line harmonic. The upper band is comparatively clean with 
no visible features. 



containing non-zero data points Tobs, and the total span 



of the observation (including data gaps) Tgpan — and 
{Sh{T)) is the noise spectral density averaged over the 
frequency search bands and the data set. Also note that 
Tobs is a function of both Tstart and Tspan- The two func- 
tions A{T) and B{T) are the integrals of the amplitude 
modulation factors and take into account the change of 
sensitivity of the instruments for the Sco X-1 location 
in the sky as a function of the time at which the ob- 
servation takes place (explicit expressions for A and B 
are given in H^). In our calculation of Q we take into 
account the presence of data gaps over T^span and we av- 
erage over the unknown angles i and tp- From Eq. H35|l it 
is straightforward to recognize that is simply related 
to the non-centrality parameter for a signal amplitude 
ho by: 



A{T) + BiT) Tobs 



hi 



(37) 



The parameter Q is therefore a faithful measure of the 
sensitivity of a given data set for the Sco X-1 search: 
it combines the effects of the variation with time of the 
detectors' noise level, duty cycle and the (angle aver- 
aged) sensitivity to the specified sky position. Note that 
for Tgpan ^ 1 day, the location of the source is a strong 
factor in the choice of the optimal data set; for longer 
observation times the quantity A + B becomes constant. 
By tuning the choice of the data set to exactly the Sco 
X-1 sky position we have achieved a gain in sensitivity 
^ 2 compared to selecting the data set only based on the 
noise level. 
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FIG. 8: The sensitivity of the LIGO interferometers during 
S2 for a search targeted for Sco X-1. The plots shows the 
evolution of the sensitivity quality factor Q, Eq. 1361 . as a 
function of observation start time; each point corresponds to 
a maximum observation time span of 6 hours. Due to inter- 
mittent loss of lock during the S2 run each 6 hour span can 
contain significantly less than 6 hours of data. Note that Hi's 
sensitivity appears more consistent than that of LI. This is 
due to the lower variation in the HI noise during the run and 
the 74% duty cycle compared to Li's 37%. The periodic struc- 
tures, more visible in the HI curve, are caused by the daily 
variation of the detectors' antenna pattern due to the Earth's 
spin. The crosses indicate the start times corresponding to 
the best sensitivity and therefore chosen for the analysis. 



We compute the figure of merit Q for all possible 
choices of data segments with Tgpan < 6 hours; the val- 
ues of Q for the whole S2 arc shown in Fig. |S1 The data 
sets that we select for the analysis begin at GPS time 
732059760, end at GPS time 732081371, and span 21611 
s with 359 60-second SFTs for HI, and begin at GPS 
time 730849195, end at GPS time 730867950, with 190 
SFTs and 18755 s for LI. Notice that T^pan is different for 
the two detectors (which has an impact on the choice of 
the orbital template banks for the two instruments) , and 
that the LI and HI data sets are not coincident in time; 
based on the relatively short observation time and there- 
fore coarse frequency resolution of the search the maxi- 
mum spin up/down of the source due to accretion would 
change the signal frequency by only '^0.1 frequency bins, 
which is negligible. 



D. Template banks 

We employ different schemes for placing search tem- 
plates for the isolated pulsar search and for the Sco X-1 
search. The optimal strategy for laying a filter bank is 
through a metric approach, and this is used for the or- 
bital parameter grid employed in the Sco X-1 analysis. 
By contrast, the search for signals from isolated objects 
uses a sub-optimal grid (in terms of computational effi- 



ciency, but not for the purpose of recovering signal-to- 
noise ratio) to cover the sky; a full metric approach was 
not developed for this search at the time that this anal- 
ysis was performed. 

To refer to the templates we will equivalently use the 
term "template" and "filter". 



1. Isolated neutron stars 

Two independent grids are employed: one in sky po- 
sition and one in frequency. The grid in frequency is 
uniform with a spacing A/q = 3.472 x 10~^ Hz which is 
about a factor of 8 smaller than the inverse of the ob- 
servation time. To cover the sky we choose an isotropic 
grid with equatorial spacing of 0.02 rad. Such a grid cov- 
ers the celestial sphere with just under 31500 patches of 
approximately equal surface area. The number of tem- 
plates in right ascension a at any given declination S is 
proportional to cos (5. At fixed a the spacing in S is con- 
stant, and equal to 0.02 rad. For illustration purposes 
Fig. 1^ shows an under-sampled grid of this type. 

The grid is chosen based on the maximum expected 
degradation in the detection statistic due to the mis- 
match between the actual position of a putative source 
and the template grid. This effect is measured by Monte 
Carlo simulations. 

The simulations consist of series of searches of signals 
at random locations in the sky with position templates 
uniformly randomly displaced from the signal's source 
position by between and half a grid step in both the a 
and the 6 direction. 

We base the selection of the grid size on the properties 
of the signals and the simulations are therefore performed 
in the absence of noise (e.g. see [13). The grid spacing is 
chosen in such a way that the expected loss in signal-to- 
noise ratio due to the signal-template mismatch is a few 
percent. The results are summarized in Fig.sllOlandllll 

The smaller the maximum mismatch between a sig- 
nal and a template the more correlated are the filters in 
the bank. We have estimated the effective number of in- 
dependent templates from the 2J- average loudest event 
found in single interferometer searches such as the one 
described here, in pure Gaussian and stationary noise in 
1.2 Hz sub-bands: 45.7 for LI and 41.7 for HI. These 
translate into an effective number of independent tem- 
plates which is a factor of ~ 4 and ~ 26 smaller for LI 
and HI respectively, than the actual number of templates 
in the sky grid that we are using for the search. The num- 
ber of independent templates was estimated asl/Po{J-*), 
where !F* is the measured loudest value of the detection 
statistic and Pq is the false alarm probability defined in 
Eq. H33|l . This is consistent with our Monte Carlo sim- 
ulations fFig.s [TUI and [TT|l where we can see that, with 
the same grid, 50% of the sky is covered in HI with a 
mismatch that is always smaller than 0.5% whereas in 
LI 50% of the sky is covered with a mismatch which is 
about twice as large. This means that the grid "covers" 
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FIG. 9: The sky grid that we have used is of the kind shown 
here, 25 times more dense. The angular distance between 
points at the equator in this plot is 0.1 radians. The figure is 
a projection of the grid on the Northern hemisphere. Distance 
along the radial direction is proportional to the cosine of the 
decUnation. 
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FIG. 11: Fraction of trials injected in HI where the ratio 
1.0 — J?-"obs/.7^optimai IS Smaller than the value on the a;— axis. 
Here in this plot .7-"obs is the .7-"-statistic at the point nearest to 
the signals' source location. In 99% of the trials the mismatch 
is smaller than 2%. 
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FIG. 10: Fraction of trials injected in LI where the ratio 1.0 — 
.?^Dbs/.?^optimal IS Smaller than the value on the a;— axis. Here 
in this plot J^obs is the .7^-statistic at the grid point nearest to 
the signals' source location. In 99% of the trials the mismatch 
is smaller than 4%. 



2. Sco X-1 

The analysis for Sco X-1 requires us to search over two 
orbital parameters and the gravitational wave frequency, 
as defined in a given reference frame that we choose to be 
the rest frame of the source (plus a correction due to the 
constant motion of the center of mass of the binary sys- 
tem). In order to optimally^ cover the parameter space 
we consider the metric approach introduced in js^ l in the 
context of binary inspirals and applied to pulsar searches 
in 0,133. We define the mismatch between a signal, 
described by the parameter vector A = {/q, Op, T}, and a 
template described by A -I- AA as® 



Ai(A,AA) = 1- 



7'(A,AA) 



= 5o/3WAA"AA'' + 0(AA3), 
where the power spectrum V is given by 



(38) 



P(A,AA) 



(39) 



HI data parameter space with more redundancy than it 
covers LI. 

The main reason for the difference in coverage is the 
fact that the spanned observation time of the data set 
used for the LI detector is longer than that for HI, and 
the resolution in sky position is highly dependent on the 
spanned observation time of the data set. 



^ Although we use the metric approach to lay templates in the pa- 
rameter space, we use a simple square grid which is non-optimal. 
Using a hexagonal grid would reduce the number of templates 
by ~ 30%. 

® Note that we use the power spectrum V to define the mismatch 
^ and therefore the metric gap, but we use the jF-statistic in the 
actual search. As the .7-"-statistic is an optimally weighted sum 
of power spectra we would expect that the template bank is also 
as effective used with the .^"-statistic as for the power spectrum. 
This was tested with extensive software signal injections. 
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A$(t; A, AA) = $(t; A) - $(t; A + AA) is the difference 
between the signal and template phase and a, P = 0,1,2 
label the search parameters (we follow the convention 
that the index labels frequency). The metric on the 
parameter space is given by [loj 



gc^p = (d^A<i>dp A$) ~ (9„A$) {dpA<^) 



(40) 



where 9aA$ = 9A$/aAA" and is evaluated at AA = 0, 
(...) stands for the time average over Tspan, and the 
gravitational wave phase is defined in Eq. (|29|l . Treat- 
ing the frequency as a continuous variable (we discuss 
later in this section the consequences of the fact that the 
frequency is in practice discrete) and "projecting out" 
the search frequency dimension of the metric yields a 
2-dimensional reduced metric only on the orbital param- 
eters, 



9Qj9kQ 
90Q 



(41) 



By using the metric 7jfc we take advantage of the correla- 
tions between the frequency and the orbital parameters 
— so that a mismatch in orbital parameters can be com- 
pensated by a mismatch in frequency — and we therefore 
reduce the number of orbital templates required to cover 
the parameter space. In the actual analysis, we carry 
out a coordinate transformation from a-p and T to two 
"search coordinates" in order to obtain constant spacing 
and orientation of the filters over the whole parameter 
space, which simplifies the numerical implementation of 
the grid. 

The frequency, however, is sampled discretely and 
therefore cannot compensate exactly a mismatch in or- 
bital parameters: it produces a mismatch /ij only in / 
between a template and a signal. In order to choose the 
appropriate frequency spacing A/ we consider the goo 
component of the metric gap — we treat de facto / as a 
1-dimensional uncorrelated dimension in the parameter 
space — and obtain: 



(42) 



^ span 



The optimal method for dividing up the total mismatch in 
detection statistic is to split them up equally amongst the 
dimensions of a parameter space. Therefore to achieve 
the required maximal overall mismatch of 10% we use a 
6.6% mismatch in orbital templates and a 3.3% mismatch 
in frequency. This sets the frequency resolution to A/ = 
l/(5Tspan) corresponding to 9.285051 x lO'^ Hz for HI 
and 1.066382 x lO"'^ Hz for LI. 

The grid in the 2-dimcnsional space of the orbital tem- 
plates is computed using ^jk, Eq. H41|) . In practice, 
the filter spacing is determined primarily by Tgpan and 
the grid orientation is determined by the location of the 
source in its orbit during the observation. These effects 
are shown in Fig. 1121 although the observation spans dif- 
fer by only ~ 15% between LI and HI, the density of fil- 
ters for HI is twice that for LI. as the number of orbital 
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FIG. 12: Here we show a small section 1/50*'') of the totai 
orbitai parameter space. The crosses and circles represent 
tempiate iocations used to search the LI and HI data sets, 
respectively. This particuiar template bank was generated 
for a maximum search frequency of 465 Hz. Note that the 
templates are not uniformly spaced in the (ap,f) parameter 
space, although they appear nearly so in the limited region 
shown here. 



templates in the regime Tgpan < P scales as T^p^^j^. The 
source location within the orbit differs by « 2.22 radians 
between the LI and HI observation periods and corre- 
lations between the two orbital parameters are therefore 
different between the two detectors resulting in template 
banks that are clearly non-aligned. One further step to 
optimize the search is to generate separate orbital tem- 
plate banks for each 1 Hz frequency sub-band, because 
the grid density increases as (/(J""')'^, where /q""' is the 
maximum search frequency. This approach allows an 
overall gain w 30% in computational speed in compar- 
ison to using a single template bank with a maximum 
frequency parameter = 624 Hz for the whole analy- 
sis. 

The number of orbital templates used for each 1 Hz 
sub-band ranges from 3391 to 3688 in the 464-484 Hz 
band and from 5738 to 6107 in the 604-624 Hz band 
for the LI analysis. The number of frequency filters per 
1 Hz band is 93,775; therefore the number of trials used 
to cover the parameter space is in the range 3.2 x lO'^ — 
5.7x10*. The corresponding numbers for the HI analysis 
are 6681-7236, 11309-12032, and 108,055, respectively, 
corresponding to a total number of trials for each 1 Hz 
sub-band in the range 7.2 x 10* - 1.3 x 10^. 



E. The single detector search 

As described in Sec. IV Bl and Fig.0]the data from each 
detector is searched over the entire parameter space, by 
computing the JF-statistic for each template in frequency 
and either position in the sky (for the all sky search) 
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or orbital parameters (for the Sco X-1 search). In both 
cases we store resuhs only for those templates that yield a 
value of the detection statistic that exceeds the threshold 
2.?^thr = 20. This choice is based on limitations on the 
size of the output files of the search. 

Our template banks are highly correlated; thus in or- 
der to decrease the number of frequency templates that 
we store, we treat as correlated the template frequencies 
which are sufficiently "close to each other". In particu- 
lar, we do not register as separate templates, templates 
which differ only by such frequencies. The frequency in- 
terval that defines how close frequencies have to be in 
order to be ascribed to the same template is estimated 
based on the full width at half maximum of the J^{fo) 
curve for a representative sample of the parameter space 
and for random mismatches between signal and template, 
as would occur in an actual search in the presence of a 
signal. The resulting frequency intervals are a few times 
10"'' Hz. 

The following information is stored for each template 
above threshold: the frequency /o at which the value of 
J- is maximum, the values of a and 6 for the template, 
the total width in search frequency bins of the points 
associated with the maximum, the mean value and the 
standard deviation of J- over all those points and the 
value of 2J^ at the maximum. The same information is 
stored in the case of the Sco X-1 search, with the orbital 
parameters Op and T instead of the sky position param- 
eters a and i5. 

The computational load for the searches is divided 
among independent machines, each searching a small 
frequency region over the entire parameter space. For 
the isolated pulsar search each CPU analyzes a 60 mHz 
search band. The processing time for both data streams 
and the entire sky is typically about 6 hours on a 2 GHz 
class computer. The typical size of the output, after 
compression, from a single detector search is around 3 
MBytes. For the Sco X-1 search, typically an individ- 
ual machine searches 0.1 Hz. The equivalent run time 
on a 60 mHz search band on the entire Sco X-1 orbital 
parameter space is approximately 9 hours. Although the 
Sco X-1 orbital templates are fewer than the sky position 
templates, the two searches are comparable in computa- 
tional time because the Sco X-1 search uses a greater 
number of (shorter) SFTs. For the particular data sets 
selected for this analysis it should be noted that with 
Li's shorter spanned observation time and fewer SFTs, 
the computational load is primarily due to the HI search. 
The output from the search in total, including both de- 
tectors and all search bands, comprises « 1 GByte of 
results, corresponding to around 700 kBytes per 60 mHz 
band for the entire orbital parameter space. 

Figures E| and El show the distribution of val- 
ues of the registered templates for sub-bands in reason- 
ably clean spectral regions in both instruments (around 
247.1 Hz and 619.5 Hz respectively) in the top two plots, 
and in less clean regions in the HI data (around 329.6 Hz 
and 465.5 Hz respectively) in the bottom two plots. In 
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FIG. 13: The circles show the distribution of 2jr values for the 
templates registered after the single detector all-sky search. 
The solid line shows the expected distribution for Gaussian 
stationary white noise. The top two plots refer to the band 
247.06-247.12 Hz. The bottom two plots show the same dis- 
tributions for the 329.56-329.62 Hz band. The expected dis- 
tribution is dominated for high values of by an exponential 
term, as is evident from the linear behavior on a semilog scale. 
In the clean 247 Hz band, the theoretical and the experimen- 
tal distributions agree very well. 



the top two plots, the distributions of values closely 
follow the expected po{2J^) distribution, Eq. (|^ . This is 
not surprising in regions free of evident disturbances, as 
already shown in [13 . Note that the highest 2JF values in 
the clean bands (top two plots) are higher in Fig. 1 141 than 
in Fig. El This is due to fact that the Sco X-1 search 
has more templates than the all-sky isolated search. 



F. Coincidence analysis 

The next stage of the analysis compares the two lists 
of values of 2!F that lie above the threshold 2J^tiir = 20 
compiled for each detector. We require that given a tem- 
plate in LI, say, there exists a template in HI such that 
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FIG. 14: The circles show the distributions of 2T values 
derived from the single detector Sco X-1 searches. The top 
two plots show the distributions for a clean sub-band in both 
detectors, 619.0-620.0 Hz. The solid curve represents the the- 
oretical expected distribution. The two bottom plots show the 
same distributions for the sub-band 465-466 Hz. In this band 
the HI results are dominated by large values of 2F . 



their locations in parameter space are consistent with a 
physical signal having triggered them both. If this is the 
case, the relevant values of the detection statistic are kept 
(the two filters are regarded as "in coincidence"), other- 
wise they are "rejected" and removed from the lists. This 
procedure is identical for both searches, but the consis- 
tency criteria are different due to the different signals 
that are searched for. This strategy is effective at re- 
ducing the false alarm rate if the noise in the two data 
streams is uncorrclated. In practice, the data arc also 
populated by a forest of lines present both in LI and HI, 
such as 16 Hz harmonics from the data acquisition system 
and the 60 Hz power line harmonics, and this procedure 
does not eliminate them. However it does eliminate the 
non-Gaussian uncorrclated outliers which also are in the 
data. We find that the typical sensitivity improvement 
in ft-o resulting from the coincidence stage is comparable 
for both searches and in the range 10% - 20%, depending 



on the frequency sub-band. 

An additional criterion to identify coincident templates 
could be based on comparing the values of IT produced 
by the two filters; however, as 2T is already maximized 
over the nuisance parameters "0 and t, and the integra- 
tion time of the analyses is shorter than 1 day, it is in 
practice difficult to introduce an "amplitude consistency 
cut" that is simultaneously stringent and safe. For this 
reason we have not included this requirement in the co- 
incidence stage of this search (sec however the discussion 
in Section IVl A|) . 

The coincident templates are then sorted in order of 
descending ioint significance. If we indicate with IThx 
and 2Th\ the values of the detection statistic for a pair of 
templates in coincidence, we define their joint significance 
as: 

s{2Tli.2Thx) = l-Po(2.fLi)Po(2^m), (43) 

where Po(2^), defined in Eq. ((SSI, is the single detec- 
tor false alarm probability for 2J-", under the assump- 
tion that the noise is Gaussian and stationary. We con- 
sider the loudest coincident template pair as that yield- 
ing the largest value of joint significance. In practice, in 
the numerical implementation we rank events according 

to -{log[Po(^Ll)]+log[Po(-^Hl)]} with log[Po(^)] = 

log(l+^)-^. 

In the remainder of the section we provide details 
about the specific implementation of the coincidence 
stage for the two analyses. 



1. Isolated neutron stars 

The candidate events that survive the coincidence 
stage are those present in both detectors' data sets and 
lie in locations of the parameter space that are consistent 
with a common signal. For the isolated search the coin- 
cidence windows are 1 mHz in frequency /o and 0.028 
rad angular distance in position on the celestial sphere. 
These coincidence window values were derived from the 
results of the Monte Carlo simulations described in Sec. 
IV El More specifically 0.028 rad represents a mismatch 
between sky positions of at most 1 grid point. The value 
of 1 mHz is derived from the results of Monte Carlo sim- 
ulations by requiring a null false dismissal rate. 

As described in Sec. IV Dl the all-sky isolated search, 
if performed on Gaussian white stationary noise, would 
yield single- interferometer loudest templates in 1.2 Hz 
sub-bands with mean 2T values of 45.7 for LI and 41.7 
for HI. The difference in these mean values is due to the 
different data sets used for the two searches (the time 
spanned by the LI data set is longer than that of the HI 
data set) and by the different location of the detectors 
on Earth and to the non-uniform antenna pattern of the 
detectors. In this search, after having excluded outliers 
with 2jr > 100, we measure mean values of the loudest 
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FIG. 15: These plots show loudest 2T values in the single 
detector searches in each 1.2 Hz sub-band. 



templates of 52.2 for LI and 46.6 for Hl"^. This corre- 
sponds to an increased level of spectral contamination in 
the real data with respect to Gaussian stationary noise. 
This is not surprising at all — even a simple visual inspec- 
tion of the spectra reveals that they are contaminated by 
several "lines" (also see the discussion in Sec. VI D in 

[13). 

After the coincidence step the mean value of IT for 
the loudest event is 39.5 for the LI data and 32.2 for HI. 
If one compares these values with the mean values before 
coincidence, 52.2 for LI and 46.6 for HI, one recognizes 
that the coincidence step yields an improvement in ho 
sensitivity of 15% and 20% for LI and HI respectively 
(remember from Eq. H35|l that 2!F oc h^). In Gaussian 
stationary noise the expected improvement is 11% and 
17% for LI and HI respectively. Thus, and again not 
surprisingly, the coincidence step plays a greater role on 
real data, which is affected by uncorrelated non-Gaussian 
disturbances. 

Figures El and El show the values of the detection 
statistic for the loudest events before and after coinci- 
dence, respectively. Figure ITTl shows the distributions of 
2!F for the loudest coincident templates. FigurelTHlshows 
the distribution of loudest coincident templates over the 
entire sky for all the 1.2 Hz sub-bands. A higher concen- 
tration of templates is apparent at the poles. This is to 
be expected since the poles are the regions from where a 
monochromatic signal would be received by our detectors 
at a nearly constant frequency. In other words, spectral 
artifacts at fixed frequency are consistent with sources 
close to the poles, during our observation time. 



^ The > 100 cut has been made only when computing the 
mean values reported above in order to eliminate large outliers 
that would have dominated the mean; see Fig.s lTOl and ll?! 



FIG. 16: These plots show two different views of the values of 
the detection statistic 2T of the loudest coincident template- 
couples (one for every 1.2 Hz sub-band) from the isolated 
pulsar search. The mean value between the threshold and 
2^ = 100 is 39.5 for the LI data and 32.2 for the HI data. 



L1 



E 

5 30 

° 20- 

<n 

± 10- 



nnnnn„nnnnnn„, 



100 

2F 



E 

Q , 




HI 



2F 



FIG. 17: These plots show the distributions of the values 
plotted in Fig. [TH The mean value of these distributions de- 
pends on the volume of the parameter space that the search 
extends over. In this case it is the whole sky in 1.2 Hz fre- 
quency sub-bands. 



2. Sco X-1 

In the Sco X-1 analysis we identify coincident tem- 
plates in LI and HI by using the metric on the rele- 
vant parameter space; see Section IV D 21 A real signal, 
if present and of sufficient amplitude, will trigger tem- 
plates in both detectors. These templates will be in close 
proximity in parameter space but not necessarily identi- 
cal. The procedure starts by taking each template above 
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FIG. 18: These figures show the location in the sky of the 
loudest coincident template in each 1.2 Hz sub-band. The 
circles show the templates found with the LI search and the 
diamonds show the coincident templates found with the HI 
search. 



threshold from the LI detector and first testing for or- 
bital parameter consistency with templates from the HI 
detector. We do this using a property inherent to the or- 
bital template bank, which is schematically illustrated in 
Fig. El We know that a signal will return at least 90% of 
its optimal detection statistic when processed using the 
"closest" template (in the absence of noise), where the 
distance between templates is defined by the metric. We 
can therefore identify a rectangular region around the LI 
filter in the 2-dimensional plane of orbital parameters. 
This region is easily calculated from the metric used to 
place the templates, and we would expect the true signal 
parameters to lie within it. We can now repeat this pro- 
cess for each HI template and construct, based on the 
metric associated with HI, the region covered by the HI 
filters (i.e. the filters that are associated with a value 
of 2T above threshold). The test now becomes a sim- 
ple matter of checking for any overlap between the re- 
gion covered by the LI filter under scrutiny and the HI 
filters. Overlap implies a possible consistent signal loca- 
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FIG. 19: Here we represent a small region of orbital param- 
eter space. The crosses represent orbital templates used in 
the LI search and the circles represent those used in the HI 
search. Note that the template orientation and spacing are 
not the same in each template bank. The bold x represents 
the location of a candidate event in the LI detector, and the 
dashed rectangle surrounding it represents the area within 
which a signal must lie if this template is the closest to the 
signal. The filled circles represent those templates in the HI 
detector that are possible coincident candidates. They are 
identified due to the overlap between their respective dashed 
rectangular regions and the LI candidate event dashed region. 
In the S2 search there are on average 12 orbital templates in 
HI that are "consistent" with each template in LI. 



tion (in orbital parameters) capable of triggering a filter 
in LI and HI from the two candidate events from a single 
common signal. In this analysis there are on average 12 
orbital templates in HI that are "consistent" with each 
template in LI. The process that we have just described 
is then repeated for all the LI filters. 

So far we have considered only the orbital parame- 
ters. The second stage to identify filters in "coincidence" 
is to test for frequency consistency amongst those fil- 
ters that have survived the previous test. The use of 
the projected metric described in Sec. IV D 21 exploits the 
correlations between the orbital templates and the grav- 
itational wave frequency to reduce the overall number of 
filters. Doing so allows greater differences between the 
true and detected source orbital parameters and greater 
differences between the true and detected gravitational 
wave frequency. Using Monte Carlo simulations it has 
been possible to define the maximum separation between 
a signal's true and detected frequency. This separation is 
very strongly governed by the spanned observation time. 
For the data sets chosen the maximum separation for 
LI was found to be 2.158 x 10"'' Hz and for HI was 
1.773 X 10"'* Hz. This corresponds to a maximum sep- 
aration of 3.931 X 10"^ Hz between a candidate in LI 
and HI in order to be consistent with a common signal. 
This is the frequency coincidence window that we have 
chosen. Note that it is equivalent to ±40 frequency bins 
in the HI search. 

If a pair of candidate events is found to be consistent in 
both orbital parameter space and frequency space then 
they are classed as a coincident event. Note that a single 
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FIG. 20: Here we show the effect that the coincidence anaiysis 
has on the loudest measured detection statistic in each 1 Hz 
sub-band within the Sco X-1 parameter space. The solid black 
curves represent the loudest coincident 2J^ values. The dashed 
gray curves represent the values before the coincidence 
analysis. Note that in clean sub-bands there is a reduction of 
~ 1.4 in the loudest value. 



candidate event in the LI detector can have many possi- 
ble coincident pairs in the HI detector (and vice versa). 

The power of the coincidence analysis is shown in 
Fig. 1201 where the effect of the coincidence constraint is 
seen to reduce the values of our loudest events. Before 
coincidence the average value of for the loudest events 
(excluding three 1 Hz sub-bands that contain major spec- 
tral disturbances: 465-466 Hz and 479-481 Hz) was 40.8 
for LI and 45.4 for HI. After coincidence this becomes 
28.6 for LI and 33.5 for HI which corresponds to an im- 
provement in ho sensitivity of ~ 16%. This is broadly 
consistent with the results obtained for the isolated neu- 
tron star analysis. In Fig. |^ we show the location of 
the coincident templates in the orbital parameter plane 
that produce the loudest event in each of the 40 1-Hz 
sub-bands. 



FIG. 21: The locations in orbital parameter space of the loud- 
est events found in each 1 Hz sub-band. The crosses represent 
events found in LI and the circles represent events found in 
HI. 



G. Upper limits 

In every ~ 1 Hz sub-band an upper limit on the ampli- 
tude of the gravitational wave signal from either a popu- 
lation of isolated neutron stars or from Sco X-1 is placed, 
based on the loudest coincident event found in that band 
during the search. The procedure employed is conceptu- 
ally identical to the one used in |l5l| to set an upper limit 
on the emission from J1939+2134, given the measured 
values of the J- statistic for that targeted search. In this 
section we describe the Monte Carlo procedure in detail. 

Let s*(/o) indicate the measured value of the joint sig- 
nificance of the loudest coincident event in the sub-band 
beginning at frequency /q. For every sub-band a set of 

injections of fake signals in the real data is performed 
at fixed amplitude ho- Each injection is searched for in 
the data and, if detected as a coincident event, its sig- 
nificance is computed. A confidence C(/io) is assigned to 
this set of injections 



Ciho) - n{ho)/N 



(44) 



with n{hf)) being the number of trials out of N in which 
the measured joint significance of the injected signal is 
greater than or equal to s* . Eq. 1)44(1 defines the upper 
limit value as a function of the confidence C . 

For every injection in a set at fixed the remaining 
parameters are chosen randomly from within the bound- 
aries of our parameter space. These include the orbital 
parameters (orbital semi-major axis and orbital phase 
reference time) for the Sco X-1 search, the sky position 
for the all-sky search, and the frequency /o and the nui- 
sance parameters ?/;, cost, and (/>o for both searches. Uni- 
form distributions are used for /o in the sub-band, V 
between and 27r, cos/, between —1 and 1, and (p^ be- 
tween and 27r. For the all-sky search the population 
of injected signals is uniformly distributed on the celes- 
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tial sphere, that is to say that a is uniformly distributed 
between and 27r and cos S uniformly distributed be- 
tween — 1 and 1, with 6 between —tt/2 and tt/2. For 
the Sco X-1 search the signal population is uniformly 
distributed across the 2 dimensional orbital parameter 
space (the parameters of which are given in Table QJ. 
The semi-major axis is selected from within the range 
1.26 to 1.62 seconds and the orbital phase reference time 
is selected from within the GPS time range of 731163028 
and 731163626. The sky position and orbital period are 
held fixed at values corresponding to the center of their 
respective ranges. For each set of injections the orbital 
eccentricity is held fixed at one of the following discrete 
values: e = 0.0, 10''', 5 x lO"*, lO'^, 5 x lO'^. 

A search over the entire parameter space is not per- 
formed to search for every injection — it is computation- 
ally prohibitive. Rather, the detection statistic is com- 
puted at the nearest template grid point with respect to 
the injected signal parameters. The nearest template is 
chosen consistently with the criteria used for laying out 
the template bank. For the Sco X-1 search the closest 
grid point is defined by the metric governing the orbital 
parameter space. For the isolated pulsar search a Eu- 
clidean measure is used. In the actual search noise might 
conspire to produce a higher value of the detection statis- 
tic at a template grid point which is not the nearest to the 
actual signal's parameters. This means that our Monte 
Carlo may slightly underestimate the detection efficiency 
of the actual search, leading to an over-conservative (thus 
still correct) upper limit. However, since our template 
bank has been chosen so that at most a few percent of 
the detection statistic may not be recovered at the near- 
est grid point due to signal-template mismatch, we do 
not expect this effect to be severe. Furthermore, a detec- 
tion/coincidence scheme based on the global properties 
of the detection statistic [s^ (far from the signal's true 
parameters) remains to be understood. 

A set of injections at fixed ho comprises at least 6000 
trials for the isolated pulsar search and 5000 trials for the 
Sco X-1 search. 

To determine the number of injections, several sets of 
10000 isolated pulsar injections have been analyzed. The 
injections were performed at a fixed strain (ho = 1.2 x 
10'^^) in a small band around 409 Hz, with sky locations 
and nuisance parameters distributed as described above. 
Figure shows the results of this analysis. We plot 
the standard deviation on the confidence as a function of 
the number of injections for seven sets of injections and 
compare it with the expected values. The plots show 
that above 5000 injections the standard deviation on the 
confidence is less than 0.2% and in agreement with the 
expectations even for small total number of injections. 

For the isolated pulsar search the following approach 
has been adopted to estimate the uncertainties related 
to a frequentist upper limit based on signal injections. 
The finite sample size of the population of signals that 
we construct by a Monte Carlo method results in fluctu- 
ations in the value of the confidence C which we measure 
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FIG. 22: A total of seven sets of ten thousand injections 
was performed in a small band around 409 Hz at a strain 
ho = 1.2 x 10~^^ . For each set the confidence was estimated 
for different numbers of injections. The plot shows the stan- 
dard deviation based on the seven estimates as a function of 
the number of injections. The dashed line shows the expected 
value of the standard deviation based on the binomial distri- 
bution with a single-trial probability of 95.6, which corre- 
sponds to the measured mean confidence at 10000 injections. 
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FIG. 23: Distribution of /j|^"j<='^'<='^ versus confidence values for 
sets of 6400 trials. The maximum spread in the confidence 
around the target 95% value is about 2% and the maximum 
spread in ho is about 10%. 360 points are plotted in this fig- 
ure, corresponding to 360 sets, each comprising several thou- 
sands injections in the 247.6-248.8 Hz sub-band. The linear 
fit from these points is C = 3.3 x 10^^ h^'^""^'"^ + 0.72. For 
C=0.95 the estimated h^^""""^ is 6.9697 x 10'^^ This is con- 
sistent with what is found in the actual analysis with a much 
smaller (order 10) set of injections: (6.999 ± 0.04) x 10"^^ 
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at fixed gravitational wave amplitude Figurcl^ 
shows a distribution of (/i™''°'^*°'^,C) values for various sets 
of injections around the target 95% confidence value in 
a reasonably clean sub-band of the data. Close to the 
target confidence the relationship C {h'^^'^'^^'^'^) is well de- 
scribed by a linear relationship. In order to estimate this 
relation we perform between 6 and 15 sets of injections. 
Each set is composed of at least 3200 injections and yields 
a value for (/iq"^'''^*'''*, C). The linear relation is then es- 
timated from these (/ig'^^°'^*°"^, C) points with a standard 
best fit technique. We define ft,Q^^° as the value of 
yielding C = 95% according to the fitted linear relation. 
From the fit we estimate the ilcr (ho, C) curves and from 
the intercepts of these with C = 95% the uncertainties 
on ft.Q^^°, which we expect to be a few percent. 

In the Sco X-1 search each set of 5000 injections is 
divided into 10 subsets, each containing 500 injections. 
The confidence Ci(/i™''°'^*°'^) is calculated for each subset 
of injections where the index i = 1, . . . , 10 labels each 
subset. Values of ho,i{C = 95%) are obtained by inter- 
polation between the two values of Ci{h™^°'^^° ) closest 
to 95% within a given subset. The final value of hl^'^" 
is calculated as the mean of /io,i(C = 95%) and the un- 
certainty in this quantity is taken as the standard error 
in the mean, cT{hf''^'') = ^(t2/(10- 1), with cr^ the vari- 
ance of the hQ^i{C = 95%) sample. This approach has 
typically yielded uncertainties in the values of the upper 
limit of - 1 - 3%. 

An additional and larger uncertainty arises from the 
instrument calibration, which varies with time and de- 
pends on the detector and the frequency band. Over the 
entire parameter space, we estimate that for the data 
sets used in the analyses presented here, the uncertain- 
ties amount to 11% and 9% for the isolated neutron star 
and Sco X-1 analysis, respectively. These estimates are 
conservative. 

The tables in [s^ detail all the upper limit results. 
The uncertainties associated with the upper limit Monte 
Carlo procedure are reported separately from the cali- 
bration uncertainties and typically they are smaller than 
the latter. 



H. Validation: Hardware injections 

Signals can be injected into the instrument via the ac- 
tuator, by physically moving the mirrors of a Fabry-Perot 
cavity to mimic a gravitational wave signal. Hardware in- 
jections are designed to give an end-to-end validation of 
the data analysis pipeline, including some, but not all, 
components of the calibration. Toward the end of the S2 
run, two simulated isolated pulsar signals were injected 
into the data. We denote the two pulsars PI and P2 and 
give their parameters in Table HTl 

The data sets were prepared using the final S2 cali- 
bration version and consist of 17 30- minute SFTs 
for HI and 14 30-minute SFTs for LI. We performed a 



PI 



P2 



fo (Hz) 1279.123457 1288.901235 

/ (Hz/s) -10"** 

a (rad) 5.147162 2.345679 

S (rad) 0.376696 1.234568 

ip (rad) 

cos L 

$0 

To (sec) 733967667.12611231 733967751.52249038 

ho 2 X 10"^^ 2 X 10"^^ 

SNR HI (exp./meas.) 17/18 34/35 

SNR LI (exp./meas.) 20/22 21/22 

TABLE II: Parameters of the two hardware injected pulsars. 
See Eqs. (1251 and 12611 for the definition of the parameters. 
To is the reference time for the initial phase, in GPS seconds 
in the SSB frame. 



targeted search to look for the pulsar signals. 

The expected SNRs are computed using a noise esti- 
mation technique that accounts for the amplitude mod- 
ulation of the signal throughout the observation time. 
The results are shown in Table ITll The agreement is very 
good. The measured SNRs, however, are systematically 
somewhat larger than expected. This is probably due to 
a small systematic error in the calibration. The differ- 
ences between the expected SNR values shown here, and 
those quoted in |l6j . arise primarily from differences in 
the lengths of observation times used to make the esti- 
mate. In |l6l |. a nominal observation time of 12 hours 
was used. This is the length of time during which the 
hardware injections were performed. Here we have used 
the actual science observation time which is shorter, re- 
flecting science quality data and calibration quality flags 
based on which we discard data as not reliable enough to 
be included in an astrophysical search. 



VI. RESULTS 

In this section we present the results of the analysis 
performed using the pipeline shown in Fig. ^ and de- 
scribed in Sectional We first discuss the results regard- 
ing the all-sky search for isolated neutron stars and then 
turn to upper limits on radiation from Sco X-1. 



A. Isolated neutron stars 

Figure 1^ shows the 95% upper limits on ho for ev- 
ery 1.2 Hz wide sub-band over the whole sky. The val- 
ues of the frequency refer to the lower extremum of each 
sub-band. (The upper limit values, along with their esti- 
mated uncertainties, may also be found in tabular form 
in \5(^.) The circles around the upper limit dots mark 
points in the 90th percentile in joint significance. About 
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FIG. 24: Upper limits based on the loudest template over the whole sky in 1.2 Hz sub-bands. The vertical stripes mark the 
sub-bands containing known spectral disturbances. The circles mark the 90th percentile most significant results. The squares 
indicate that the values of the detection statistic in the two detectors are not consistent with what one would expect from an 
astrophysical signal. 



2/3 of these points are also in the 90th percentile for 

1,95% 

"o 

About one quarter of the 90th percentile significance 
points lie in sub-bands influenced by spectral distur- 
bances (points with circles on shaded bands in Fig. I24|l . 
Most of the remaining points can be immediately at- 
tributed to non-astrophysical sources because the ratio 
of the JF-statistic values in the two detectors is cither too 
large or too small to be consistent with being due to the 
same signal. These points are indicated in Fig. [21 by a 
square. There remain 6 points which are in the 90th per- 
centile in significance and cannot be excluded based on 
the ratios of the JT-statistic values. They appear at the 
frequencies 160.0 Hz, 466.79 Hz, 546.03 Hz, 564.02 Hz, 
626.80 Hz, and 700.51 Hz. The 160.0 Hz frequency coin- 
cides with the 10-th harmonic of 16.0 Hz, a key operating 
frequency of the data acquisition system. We are thus 
confident that the origin of this outlier is instrumental. 
The points at 466.79 Hz, 626.80 Hz, and 700.51 Hz are 
due to lines only in the LI instrument which have dis- 
appeared in science runs subsequent to S2. This check 
suggests that the outliers are of instrumental origin. The 
546.03 Hz and 564.02 Hz points are due to lines which 
clearly appear only in HI. However the lines are present 



in the S2 nm and in later science runs. The amplitude 
of both lines decreases with increasing sensitivity of the 
instrument, dropping by a factor of 10 (in noise power) 
as the sensitivity increases by a factor greater than 5. 
This indicates a behavior which is not consistent with 
the model of the signal that we are searching for here 
and suggests that these lines are of instrumental origin. 
Figure shows the average power spectral density in 
both detectors in the frequency regions where these out- 
liers are located. 

Unlike what is described in |l3l | no frequency band is 
excluded from the upper limit analysis due to it being 
contaminated by known noise artifacts. This results in 
extremely loud events in some sub-bands: those contain- 
ing the 60 Hz power line harmonics, the LI calibration 
line (at 167 Hz), the violin modes of the suspension wires 
of the test mass (in the 340-350 Hz region) and the var- 
ious oscillator harmonics at multiples of 36 Hz together 
with the beating of the 0.74 Hz pendulum mode of a test 
mass against the oscillator line (in the 220-335 Hz re- 
gion). In the case of the 179.4 Hz sub-band containing 
the 180 Hz power line harmonic, the spectral disturbance 
is so strong that the upper limit Monte Carlo does not 
converge to an upper limit h^^'^" value. The sub-bands 
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FIG. 25: The average power spectral density in the two detectors in the frequency regions around the five not immediately 
explained outliers from the search for isolated pulsars. The width of these disturbances is sufficiently small that they could not 
be discarded as of non-astrophysical origin based on this. A peak in power spectral density is clearly visible only in one of the 
detectors, but the measured values of the detection statistics are not inconsistent with an astrophysical signal, albeit a rare 
one. As explained in Section [Vl Al we very strongly suspect these excesses of power are not due to a continuous wave source 
because of their inconsistent amplitudes in science runs subsequent to S2. 



marked by shaded vertical stripes indicate frequencies 
where known spectral artifacts are present. 

The upper limit values presented here are in broad 
agreement with what is expected and we consider this a 
further validation of the analysis pipeline. We have run 
the pipeline presented here on Gaussian stationary noise 
and empirically derived a formula for the expected /ig^'^" 
as a function of the noise level in the detectors: 



with (Sh) being the average noise level over 1.2 Hz and 
over the observation time in every detector and then aver- 
aged over the detectors: this quantity is shown in Fig. 
We would like to stress that Eq. (|45|l refers to this partic- 
ular analysis and pipeline. The expected upper limits for 
Gaussian stationary noise arc plotted against the mea- 
sured ones in Fig. [211 It is clear that in regions where 
the data is not Gaussian and stationary Eq. H45|l does not 
predict correctly the values that wc measure and the dis- 
crepancy between the prediction and the measured value 
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FIG. 26: Top plot: the circles are the upper limits based 
on the loudest template over the whole sky in 1.2 Hz sub- 
bands. The solid line is the expected upper limit under the 
assumption that the noise in the detectors is Gaussian and 
stationary. Bottom plot: ratio of the measured upper limit 
values to the expected upper limit values. The circles indicate 
the 90th percentile values of the measured to expected upper 
limit ratio. The crosses mark the upper limit values for which 
the ratio of the detection statistic values in the two detectors 
is not consistent with what we would expect from a signal. 
In both plots the shaded regions indicate frequency bands 
affected by known spectral disturbances. 



depends on the details of the spectral disturbance and of 
the method used for estimating the noise. This is particu- 
larly evident close to spectral disturbances, where clearly 
the noise is not white Gaussian and often not stationary 
and the predictions can even be larger than the actually 
measured upper limit value (see the points below 1 in the 
lower plot of Fig. l26|l . However the ratio of the measured 
upper limits to the expected one never exceeds 4.4 and 
the 90th percentile level in this ratio is 1.7. 



B. Sco X-1 

The upper limits on gravitational waves from Sco X-1 
are summarized in Fig. 1271 (more details are provided in 



[5^): we show /iq^^° over 1 Hz sub-bands in the range 
464-484 Hz and 604-624 Hz, assuming that the source is 
in an exactly circular orbit. We would like to stress that 
these limits apply to a source whose orbital parameters 
lie in the region reported in Table corresponding to 
1-a errors. The typical value of h^^'^" is ~ 2 x 10"^^ 
over the whole analyzed 40 Hz band, with the exception 
of a band w 2 Hz around 480 Hz which corresponds to 
one of the strong harmonics of the 60 Hz power-line, cf. 
Fig. El In this region the upper limit is Hq^^" ?a 10~^^. 
Such values are consistent with the sensitivity estimates 
shown in Fig.|21 which were derived under the assumption 
of Gaussian and stationary noise and include a number 
of approximations to quantify the effects of each stage 
of the pipeline considered in this search. Through the 
statistical modelling of the pipeline we are able to express 
the expected /iq^^° upper limit as 



C%(/)^28 



{Skif)) 
6 hrs 



(46) 



where {Sh{f)) is the noise level over the 1 Hz sub-bands 
averaged over the observation time, the frequency band 
and the detectors. 

In Fig. I2H1 we show the value of as a function of 
C for selected frequency sub-bands. We have considered 
both "quiet" and "noisy" spectral intervals, but have re- 
stricted this analysis to only four 1 Hz frequency sub- 
bands due to computational burdens. Figure 1281 shows 
that for C = 0.99 (0.5) the upper limit on ho would be a 
factor sa 2 larger (sa 2 smaller) than h^^'^". 

So far we have restricted the discussion of the upper 
limits to the case of an exactly circular orbit. This is the 
model that we have assumed in building the templates 
used in the analysis. As we have discussed in Sec. lIVir2l 
the orbital fits of the optical data are consistent with 
6 = 0, which is in agreement with the theoretical ex- 
pectations inferred from evolutionary models. However 
present observations are not in a position to constrain 
the eccentricity to e ;S 10~*, which would introduce sys- 
tematic losses of signal-to- noise ratio smaller than w 0.1, 
the value of the mismatch adopted for this search. It 
is therefore important to explore the consequences of a 
(unlikely but possible) non-zero eccentricity of the Sco 
X-1 orbit on the results reported so far. The pipeline 
that we have developed allows us to quantify this effect 
in a fairly straightforward way: the Monte Carlo soft- 
ware injections used to set upper limits are performed 
again by drawing signals from a population of binaries 
where now the eccentricity is set to a (constant) value 
e 7^ 0. The orbital parameters and the frequency are 
chosen randomly exactly as in the case for a circular 
orbit. We detect the signals from eccentric orbits with 
the search pipeline constructed with a bank of filters for 
a perfectly circular orbital model. In this way we can 
quote consistently an upper limit on Hq for e 7^ 0. We 
repeated this procedure for selected values of the eccen- 
tricity, e = 10"'', 5 X 10"'', 10-3 and 5 x 10"^ The depen- 
dence of the value of the upper limit on ho as a function 
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FIG. 27: The upper limits on the amplitude of gravitational wave radiation from Sco X-1. The plot shows the limits on ho at 
95% confidence (solid line) as a function of frequency. We report one limit for every 1 Hz sub-band and the Sco X-1 orbit is 
assumed to be exactly circular. The shaded region shows the combined errors on h^^'^" due to the injection process (typically 
in the range ±2 — 5 x 10^^**) and instrument calibration. The latter, estimated at the level ~ 9%, dominates the uncertainties. 



of the confidence for four representative frequency sub- 
bands is shown in Fig. 1281 and the upper limits over the 
whole 40 Hz region (at fixed confidence) are summarized 
in Fig. 1291 and in [53 ■ Notice that in this case we choose 
different values of the confidence depending on the ec- 
centricity of the orbit of the putative source population 
used for the injection. This stems directly from the fact 
that the detection efficiency of the search pipeline is pro- 
gressively reduced as the model of the injected signals 
differs more and more from that of the detection tem- 
plates. In other words we suffer from systematic losses 
of signal-to-noise ratio due to the fact that the templates 
are not properly matched to the signal: for e ^ 10"'^ the 
fitting factor of a filter generated by modelling Sco X-1 
as a circular orbit binary is < 0.9. Indeed, regardless of 
the strength of the injected signals, the pipeline is unable 
to detect at least 95% of them, see Fig. We find that 
for e = 10""^ and 5 x 10""^ the pipeline has a maximum 
detection efficiency in the ranges ~ 90-95% and ~ 50- 
70%, respectively. As a consequence, for e — 10^^ we 
report /ig^'^", and for e = 5 x 10~^ wc consider /iq^*^", be- 
cause across each of the 1 Hz sub-bands we have achieved 
at worst an 88% and 50% confidence, respectively (see 
Fig. 1231 ■ On the other hand, for e < 5 x 10""^ the sys- 
tematic loss of signal-to-noise ratio is small or even neg- 
ligible and we can quote upper limits at 95% confidence. 
Wc find that, as expected, the values of /ig^'^" for e = 
and e = 10"^ are essentially identical. For e = 5 x 10""*, 



/ig^^° is about 50% bigger than in the case e = 0. 

Considering the limited sensitivity of the present anal- 
ysis (see Fig.s O El ^nd I29|) with respect to the astro- 
physical predictions, Eq. (|24|1 . we have not followed up 
{e.g. using a longer integration time Tspan) regions of the 
parameter space that yielded particularly large values of 
J-. Such a follow up would be computationally very in- 
tensive and the fact that we are targeting a continuous 
gravitational wave emitter allows to go back to the same 
parameter space in the future, exploiting higher sensitiv- 
ity, better quality data, and a more sensitive search al- 
gorithm. This work is already in progress. It is however 
important to establish that the results that we have ob- 
tained do not show any obvious unexplained feature and 
are qualitatively consistent with the expectation that no 
signal is present in the data set at the sensitivity level 
of the search. In Fig. (201 we show the distribution of the 
parameters that characterize the filters in coincidence in 
LI and HI for a representative frequency band, 614 - 
614.2 Hz. Due to the high correlation of the templates 
used in the analysis one would expect a cluster in pa- 
rameter space of filters in coincidence, were a real sig- 
nal present. Considering the coincident filters in the 3- 
dimensional search space (/q , Cp , T) and projecting them 
onto the plane (op, T), no particular structure is evident, 
with coincident templates evenly distributed across the 
plane. This is also broadly consistent with the distribu- 
tions of coincident templates that we have obtained by 



30 











































464^65 Hz| 



o 0.; 













































: 4BM70Hz| 



2 X 10"" 



4x10"" 10"" 


































■ 614-615 Hz| 



2 X 10-=^ 10-== 



FIG. 28: Examples of how the confidence in the upper hmit 
on ho from Sco X-1 scales with the amplitude of injected sig- 
nals. Here we show 4 plots each corresponding to a different 
1 Hz sub-band. The sub-bands 464-465 Hz, 469-470 Hz, and 
614-615 Hz show the confidence across the range of typical 
injected ho values (10~^^ - 4 x 10"^^). The sub-band 479- 
480 Hz contains a large spectral disturbance and the injections 
are appropriately increased in amplitude in order to achieve 
the required confidence. The 5 curves in each plot represent 
the confidence associated with different values of orbital ec- 
centricity. The circles represent e = 0, the plusses represent 
e = 10"*, the crosses represent e = 5 x 10~*, the squares 
represent e = 10""^ and the triangles represent e = 5 x 10~^. 
Note that for the 6 hour observation, the pipeline is as sen- 
sitive to signals with values of e < lO"'' as it is to circular 
orbits. For signals with e — 10~^ the pipeline has a maxi- 
mum detection efficiency of ~ 90-95%. For e = 5 x 10~^ 
maximum detection efficiency reaches only ~ 50-70%. 



the 



performing Monte Carlo simulations of the entire search 
pipeline on stationary and Gaussian noise. There is how- 
ever some structure in the plane (ap,T) that is deter- 
mined by narrow spectral disturbances and accounts for 
the large outliers in the values of the detection statistic 
at the end of the analysis pipeline. In order to explore 
this, it is useful to project the same 3-dimensional pa- 
rameter space onto either the (/o,ap) or (/o,r) plane. 



Structures are now clearly visible consisting of "stripes" 
of events at approximately (but not exactly) constant 
frequency. These are caused by small narrow spectral 
features present in the data that produce relatively large 
values of the ^-statistic for a number of orbital tem- 
plates. Due to the short coherent integration time, such 
disturbances arc not averaged out by the demodulation 
process and arc registered in the single detector search. 
We find that they are very common in the output of the 
single detector search; the coincidence stage of the anal- 
ysis allows those "stripes" of events that exist in both 
detectors at approximately the same frequency to sur- 
vive the entire pipeline. 

We have so far reported the results of the analysis as 
upper limits on the signal amplitude /iq. We can now 
re-cast them as upper limits on the ellipticity e of the 
neutron star, taking the distance of Sco X-1 as = 2.8 
kpc; see Tabled Using 
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(47) 

and the canonical value for the principal moment of in- 
ertia Izz = 10^^ gcm^ we obtain: 



,95% _ 



4.0 X 10-4 - 3.6 x 10-3 for e < 10"'* , (48) 



over the frequency band. The previous result can be 
generalised to the case of a more pronounced non zero 
eccentricity; e.g. we obtain 



,88% 



5.1 X 10" 



3.7 X 10- 



for 



10- 



(49) 



Despite being far from astrophysically interesting, 
Eq. (|48|l and H49|) represent the first direct measurements 
of the ellipticity of the neutron star in Sco X-1 in the rel- 
evant frequency band. 



VII. CONCLUSIONS 

We have presented here results from two coherent wide 
parameter space searches for continuous gravitational 
wave signals. A subset of data from the second science 
run of the LIGO instruments was analyzed, the data cho- 
sen to maximize the sensitivity of the search. Two dif- 
ferent astrophysical searches were performed: an all-sky 
search aimed at signals from isolated neutron stars and 
an orbital parameter search aimed at signals from the 
neutron star in the binary system Sco X-1. Both searches 
also cover a wide range of possible emission frequencies: 
a 568.8-Hz band for the isolated pulsars search and two 
20-Hz bands for the Sco X-1 search. 

The sensitivity of these analyses makes the detection 
of a signal extremely unlikely. As a consequence the 
main goal of the paper is to demonstrate an analysis 
method using real data, with a pipeline considerably 
more complex than any other coherent searches previ- 
ously performed. More importantly, this coherent search 
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FIG. 29: The upper limits on ho as a function of frequency for different orbital eccentricities for the Sco X-1 search. Notice that 
for sake of clarity, we have adopted a different scale on the vertical axis than that used in Fig. 1271 as a consequence the upper 
limits in the frequency region 479 - 481 Hz are not shown because they are off scale. Due to the systematic loss of signal-to-noise 
ratio due to signal-template mismatch for populations of signals from sources in eccentric orbits, the confidence C at which 
the upper limit h^' is computed is different depending on the value of eccentricity e; see text for a detailed justification. We 
report hf^" for e = 10"*, hl'^'^" for e = 5 x 10"'', /ig**^" for e = 10"^ and /ig"'^" for e = 5 x 10"^ The shaded areas represent 
the combined errors due to the injection process and the instrument calibration. 



will be deployed in a hierarchical analysis scheme. The 
first step of a hierarchical analysis sets the ultimate sen- 
sitivity of the search: candidates that do not survive the 
first threshold are lost. It is thus crucial to employ the 
most sensitive possible technique in this first step. The 
coherent method described in this paper provides an im- 
plementation of such a first step. 

Overall, hierarchical approaches are expected to 
achieve optimal sensitivity at constrained computational 
resources. We will employ such approaches for deep 
searches on long duration and high sensitivity data such 
as those that are now being recorded by the instruments. 
With one year of data at the design sensitivity of the de- 
tectors, the improvements that we can expect by means 
of hierarchical schemes that utilize this type of coherent 
analysis as one of the steps is of order 10 with respect to 
what was presented here. 

The most constraining 95% confidence /iq upper limit 
from the all-sky search is 6.6 x 10"^^ in the band 245.2- 
246.4 Hz, reflecting the highest sensitivity of the instru- 
ments at these frequencies. This is still a factor of ~ 16 
higher than the strongest signal that we expect based on 
the optimistic (but not unreasonable) assumptions of the 
statistical argument presented in Sec. IIII B II 

The 95% confidence upper limits from the Sco X-1 
search, assuming a non-eccentric orbit, are ho ~ 2 x 
10"^^; the most stringent 95% confidence upper limits 
from the Sco X-1 search, assuming a non-eccentric or- 
bit, are ho = 1.7 x 10"^^ in the 464 - 484 Hz band and 
ho = 2.2 X 10^22 ^jjg _ 624 Hz frequency band. 
The Sco X-1 results presented here are the first direct 
gravitational wave upper limits placed on the system. 



Coherent all-sky searches for continuous signals from 
isolated stars have been performed in the past, but over 
much smaller parameter space. In an all-sky 0.76 Hz 
band search was performed around 921.38 Hz, including 
spin-down parameters in the range —2.36 x 10"^ Hz s"^ 
to -1-2.36 X 10"^ Hz s"^ Three data sets, each 48 hours 
long, were coherently analyzed and a 90% confidence up- 
per limit was placed at the level of 1.0 x 10"^^ based on 
the the cleanest of the data sets. 

In the context of a hierarchical search aimed at detect- 
ing a signal, data cleaning procedures to remove noise 
artifacts are likely to be employed. This is in contrast to 
the approach used in the analysis reported in this paper 
where no cleaning at all was considered. Here we wanted 
to investigate the effect of noisy data segments with a 
variety of artifacts on the output of a search pipeline; we 
have purposely kept the "bad" data in the presentation 
of the results for illustration purposes. This is relevant to 
future searches because the upfront "cleaning" of known 
noise artifacts does not guarantee that longer observa- 
tion times will not uncover unknown periodicities. These 
must then be cither identified as of instrumental origin, 
as done here, or followed up. 

Longer observation times mean a higher resolution in 
parameter space and higher computational costs. It is 
thus important to lay the template banks in a way that 
takes advantage of the correlations in parameter space. 
The Sco X-1 search presented here already does this. In 
future work the metric approach will be used also for 
template placement for the searches for signals from iso- 
lated sources. However, longer observation times and 
targeting different LMXB systems will require more so- 
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phisticated template placement strategies than the one 
presented here for Sco X-1. For example the parameter 
space may grow to include the orbital period, the eccen- 
tricity and the spin period derivatives. For any specific 
source the number of search parameters will be defined 
by the precision to which these source parameters have 
been measured via electromagnetic observations. 

Redundant template grids produce redundant events. 
This increases the false alarm rate and in practice re- 
duces the sensitivity of the search; in fact the thresh- 
old on signal-to-noise ratio that defines the candidates to 
follow up depends by how many follow-ups one can af- 
ford with given computational resources. It is therefore 
important for future work to develop techniques to rec- 
ognize non-independent candidates in parameter space, 
rank them and keep the information on only the most sig- 
nificant. In this paper we have made the first moves in 
this direction with the algorithm that identifies as a sin- 
gle candidate values of the T statistic which are "near" 
to each other in search frequency. The concept must be 
generalized to the multi-dimensional space of the search 
parameters, and ultimately connected to the global prop- 
erties of the detection statistic over the parameter space. 
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FIG. 30: The locations of coincident templates within a 
representative subsection of the Sco X-1 parameter space. 
These plots contain events found in the frequency band 614.0- 
614.2 Hz and are representative of a "clean" search band. We 
show only events from the HI detector for clarity (LI coinci- 
dent events lie in approximately the same locations) . We show 
three 2-dimensional projections through the 3-dimensional 
search space. The first plot shows the projected orbital semi- 
major axis, ftp, versus the orbital phase reference time, T, 
and shows no obvious structure. The second and third plots 
show ftp and T, respectively, versus /o. Here we do see struc- 
ture caused by small narrow lines present in the data that 
the demodulation process has failed to smooth out over the 
relatively short 6 hour observation time. 
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